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CONTI HJED DEVELOPMENT AND CORRELATION OF 
ANALYTICALLY BASED WEIGHT ESTIMATION 
CODES FOR WINGS AND FUSELAGES 

By Joseph Mullen, Jr. 

Nielsen Engineering & Research, Inc. 

SUMMARY 

The implementation of the changes to the program for Wing 
Aeroelastic Design (WADES) recommended in part I of this work and 
the development of a program to estimate aircraft fuselage weights 
are described. The equations derived to implement the modified 
planform description, the stiffened panel skin representation, 
the trim loads calculation, and the flutter constraint approxima- 
tion are presented. A comparison of the wing model with the ac- 
tual F-5A weight material distributions and loads is given. 

The equations and program techniques used for the estimation 
of aircraft fuselage weights are described. These equations were 
incorporated as a computer code. The weight predictions of this 
program are compared with data from the C-141. 




INTRODUCTION 


Under the sponsorship of the National Research Council from 
1972 to 1974 , a program for the aeroelastic design of simplified 
conventional and multilayered composite wings for strength and 
flutter requirements (WADES) was written in order to study tech- 
niques in structural optimization. Under Contract No. NAS2-8558* 
from NASA/Ames Research Center , Nielsen Engineering 6 Research, 
Inc. was funded to incorporate this capability into ARC'S air- 
craft synthesis program, ACSYNT. In the first phase of this 
contract detailed comparisons of the estimated weight, material 
distributions, and loads of the WADES program with those of the 
F-5A/B wing were made. As a result of that comparison certain 
program deficiencies were identified (ref. 1). In the second 
phase of this contract the changes recommended to correct those 
deficiencies were incorporated into the WADES program. 

Part of the first phase of Contract NAS2-8558 was the inte- 
gration of the wing aeroelastic design program as a module of the 
vehicle synthesis program, ACSYNT. In order to fully assess the 
potential of the advanced structural technology of this program 
on vehicle weights it was decided that both wing and body weights 
should be computed. The development of a program to predict 
fuselage weights was then undertaken as part of the second-phase 
work. A survey of current programs was made to take advantage of 
available technology. This is the final report summarizing the 
modifications to the WADES program and the derivation of the 
equations used in the prediction of the fuselage weight estimates. 

The survey of computer codes was undertaken to compare the 
technology and adaptability of available computer programs for 
the estimation of wing and fuselage weights, m summary of the 
results of that survey is given here. As a result, some of the 

* 

Technical Monitor: Dr. G. N. Vanderplaats. 
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methods employed by reference 2 have been adapted for use in this 
work. The purpose of this report is to explain the adaptation of 
those methods and the modifications are required to incorporate 
them into the current programs. Rederivation of much of the work 
in terms of geometric descriptors used by the ACSYNT program was 
required. The derived equations, the general flow of the programs 
and the comparison of the computed results with actual aircraft 
data are given here. 


3 


l kt w-zf.-Vts 


SYMBOLS 
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dU,n) 
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h 

i 

K 

AK 

M 



M 


Chord length 

Depth function describing the wing thickness over 
the planform 

Discrete force or mass of i'th component 
Altitude 

Subscript or component locati . n in an array 
Weight coefficient 
Incremental weight coefficient 
Mach number 

Bending moments about x- and y-axes, respectively 
Limit load factor 
Dynamic pressure, or shear flow 
Reaction forces at wing-body interface 


R 


r 

S 

SPAN 

t (c,n) 


t/c 

v 

v 

7 . 

w 

w/s 


x,y ,z 


x/c 

X LE, X TE 

6 1' 6 2 
y 


Root chord 

Local fuselage radius 
Planform area 
Semi-span 

Thickness function describing skin cover gage over 
the planform 

Thickness to chord ratio 

Velocity 

Shear force perpendicular to wing 
Weight 

Weight per unit area 

Basic coordinates: x-streamwise , y-spanwise, 

z-transverse 

Ratio of local x distance to chord 

x/c locations from leading and trailing edges of 
structural planform 

Leading and trailing edge sweep angles 
Poisson's ratio for material 
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SURVEY OF PROGRAMS FOR THE ESTIMATION OF WING AND FUSELAC' WEIGHTS 


A survey of weight estimating programs was undertaken in or- 
der to take advantage of prior technology and available computer 
codes in the development of a structural weight estimation for 
use with the ACSYNT program. This summary does not contain all 
the programs written to perform this task, but hopefully enough 
of the codes have been examined to formulate a representative ap- 
proach to module development. 

In the evaluation of the various programs, primary emphasis 
was placed upon their applicability to vehicles of interest to 
ACSYNT and their adaptability to vehicle synthesis. Where only 
the data base was of interest, the information was considered 
primarily for use by the WADES program. The scope of the codes 
examined was limited to those programs which predict the weights 
and structural responses of wing and fuselage structural compon- 
ents. The specific criteria employed included the following: 

(a) The level of sophistication of the program was assessed 
to determine the complexity of the structural model and the com- 
putational speed of the program. This included a judgement of 
the program's potential use in either a Level I or Level II struc- 
tural analysis or design in conjunction with the ACSYNT program. 

A Level I program in these contexts is a very fast executing pro- 
gram which generally provides a Group Weight Statement estimation 
of the component weight items based on the gross parameters of 
the aircraft. A Level II program requires longer execution times 
and generally is based on structural analysis and can provide 
weight estimation for Detail Weight Statement items. It usually 
has the analytical basis to assess new technology. 

(b) Whether primary and/or secondary structural weights are 
estimated was determined. 


6 



(c) The technology level was assessed to determine the depth 
of the structural concepts used, the advanced design techniques 
used, and the ability of the program to consider multiple design 
criteria and advanced materials. 

(d) The program's adaptability for use with the ACSYNT pro- 
gram w as determined. This included an assessment of its potential 
for direct incorporation as a module of ACSYNT, its use as a data 
source (theoretical or empirical) , its use as a stand-alone pro- 
gram to be employed vehicles under consideration. Its utility to 
ACSYNT was based initially on its adaptability to the WADES, SSAM, 
or SAD programs currently in use. 

(e) The integration effort and data required from other 
modules to use the program were estimated. 

The following programs were examined in the literature. They 
are partitional according to types of analysis used to obtain the 
weight estimates: empirical, semi-empirical, or analytical. 

Em irical Weight Estimation Programs 

The most prevalent weight-estimating techniques used in pre- 
liminary vehicle synthesis are based on empirical data. Statis- 
tically correlated equations are derived to predict the weights 
of various components. The obvious pitfall is that the equations 
are only valid within the bounds of the original data. The use 
of such techniques is limited in value when considering new con- 
cepts that may lie outside the original data. 

Three programs in this class were looked at closely. Each 
of these was being used in a vehicle synthesis program. 

(1) WAVES Program - this is currently used in the ACSYNT 
program at Ames Research Center; maintained by Alice Barlow. It 
is a Level I program for conventional transport, fighter, and 
bomber aircraft types. The group weights of combined primary and 
secondary weights are computed. Little advanced technology data 
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were used in the derivation of the empirical coefficients. Ad- 
vanced technology factors are available as user input. 

(2) WAATS Program (ref. 3) - this is a Level I program in- 
cluding component weight equations for conventional transport and 
military aircraft including some component weight estimations for 
high- temperature (X-15) aircraft. Weight estimation includes both 
primary and secondary group weights. Advanced technology used in 
correlations includes some high-speed and high- temperature air- 
craft data. Programming features include a generalized form of 
the equations with coefficients input at execution time. The pro- 
gram or equations derived therein may be easily adapted to ACSYNT. 
Component correlations demonstrated no better agreement than WAVES 
equations. This program was originally written for use with ODIN 
at LRC and is also being used with EDIN at JSC. 

(3) WTSIZ (refs. 4, 5 and 6) - this General Dynamics Level 
I program calculates the basic structural component weights of 
advanced fighter, bobmer or cargo aircraft. Weight estimation 
includes both primary and secondary group weights. No explicit 
advance technology factors were used in correlations with data. 
Programming features include an interactive mode in combination 
with performance and geometry modules for use in vehicle synthe- 
sis. Balance capability is also included. The program is pri- 
marily usee in an interactive mode. It does contain the estima- 
tion of a number of subcomponent weights not available in the 
ACSYNT . cogram. 

Other statistical weight estimating methods and programs 
exist (refs. 7-14) . Only the previous three were summarized with 
any depth, since the primary interest here is to examine semi- 
empirical and analytical approaches. Some other programs in use 
are : 

HIPERAC - program used by Naval Air Development Center to 
study high-performance aircraft. 
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VASCOMP - program used by Aeronautical Systems Branch at ARC 
for short-haul CTOL-V/STOL aircraft; developed by 
Boeing Vertol. 

HESCOMP - program used by Aeronautical Systems Branch at ARC 
to study rotary wing aircraft; developed by Boeing 
Vertol . 

GASP - program used by Aeronautical Systems Branch at ARC 
to study general aviation aircraft. 

EDIN - uses the mass- and volumetric-property programs 

SSP, VAMP , VASC , WAATS, CASPER, ESPER, and APSB. 

CASP - program used by Prototype Division of Air Force 

Flight Dynamics Laboratory for vehicle synthesis 
using WTSIZ routines. 


TABLE I. - SUMMARY OF EMPIRICAL WEIGHT PROGRAM CAPABILITIES 


Capability/Program 


WAVES 


WAATS 


WTSIZ 


LEVEL 


Weight Types Computed 
Primary 
Secondary 
Subcomponent 

Technology Level 
Conventional 
Composites 
Thermal 
Input 

Program Features 

Default Values 
Optional Input 
Generalized Equations 

Vehicle Applicability 
Transports 
Fighters 
Bombers 

General Aviation 
RPV's 

High-Speed Aircraft 

Missiles 

Cargo 

ACSYNT Compatibility 
Code 

Data Base 
Stand-alone 
















Semi-Empirical Weight Estimating Programs 

The second grouping of weight estimating methods examined is 
classified as semi-empirical. These techniques are analytically 
derived explicit or integral equations which size the structure 
or structural component to satisfy a single critical design cri- 
terion. The equations are then multiplied by a statistically ob- 
tained "non-optimum" factor to compute the actual structural 
weight. Because these equations typically try to predict weights 
of such complex structural components as wings and fuselages , 
there can be 20-30 percent error resulting from such "non-optimum" 
factors. The following semi-empirically derived weight estimating 
methods were examined: 

(1) The wing and fuselage weight estimating equations by 
Shanley and Micks are given in reference 15 and 16. The explicit 
forms of the equations could be developed into Level I weight- 
estimating relationships. A numerically integrated form could be 
developed to handle geometric discontinuities within the gtiide- 
lines of Level I. The weights computed are the primary structur- 
al weights, with the secondary weights calculated by the use of 
"non-optimum" factors or empirically based secondary component 
equations. The technology level used to derive these equations 
was derived from conventional metal construction. The analysis 
assumes certain optimal failure modes in deriving sizing criteria. 
Subsequent modification of these optimal relationships may be 
incorporated to provide limited capability for advanced technol- 
ogy. No direct program was currently available, though one could 
be assembled using the ACSYNT geometric descriptors. Minimal 
data from other modules would be required. 

(2) Wing and body weights are computed as part of the TRANSYN 
program in reference 12 for transport-type aircraft. A similar 
derivation for elliptically shaped hypersonic configurations is 
found in reference 18. The transport wing weights are computed 
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using a multi- spar box beam sized on the basis of the critical 
instability modes in reference 19. The volume of material required 
for shear and bending is computed at each station and integrated 
spanwise. The structural model, however, incorrectly models shear 
flow in the wing. Torsion is included only in the carry-through 
structure, and only maximum symmetric pull-up loading is consid- 
ered. Considerations of the detail of the structural model and 
relative computational speed suggest that this program is Level I. 
Secondary weights are computed by using a "non -optimum" factor. 
Advanced technology is limited to suitable modifications of crit- 
ical instability modes. Design algorithms are limited to one- 
dimensional searches for point designs. The program follows a 
Shanley approach and does not provide sufficient improvement to 
merit direct incorporation. 

The fuselage weight computation is an adaptation of Shanley 's 
sizing procedure. Body sizing is based on maximum bending moment 
due to symmetric pull-up or dynamic landing loading conditions. 

This is also a Level I program. Secondary structural weights of 
attachments, bulkheads, etc. are estimated by a single "non-opti- 
mum" weight factor. No advanced technology has been incorporated. 
The only improvement over Mick's development is the inclusion of 
numerical integration along the body length, and consideration of 
the landing loads. No combined loading effects or area ruling 
effects are considered. 

(3) A wing-box weight predicting method is proposed by Burt 
in references 20 and 21. His method basically follows the approach 
taken by Shanley of computing the volume of material needed for 
shear, bending, and torsional requirements. Burt's contribution 
was the consideration of several types of construction (skin- 
stringer-rib, honeycomb, corrugated core, etc.) in the weight 
equations. No program exists that could be readily adapted. It 
would be a Level I program. Only the primary structural weights 
are predicted analytically. The only advanced technology concepts 
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included in his derivations were the optional construction types. 
Most of these are mid-50's construction state-of-the-art. The 
associated data would be the only portion of the analysis with 
potential for adaptation as a structural weight estimation module 
of ACSYNT . 

Analytically Derived Weight Estimating Programs 

The third grouping of weight estimating programs examined is 
classified as analytical. These programs are generally derived 
to model the actual structural concept and analyze and design the 
structure to satisfy multiple loading conditions. Secondary struc- 
tural component weights are usually included by either empirical 
or semi-empirical relationships. Analysis models fall into two 
categories: beam models and finite-element models. Design ap- 

proaches include the use of point sizing, optimality criteria, 
and mathematical optimization and most of these programs use Level 
II category analyses and design. The following analytical weight- 
estimating programs were examined: 

(1) The preliminary wing design program developed by Harold 
Switzky (Fairchild Republic Co.) is outlined in reference 22. 

This program is directed at the design of high-aspect-ratio box 
structures and uses a beam station analysis. Loads are generated 
by a vortex-lattice aerodynamic technique and scaled during the 
design to account for the changes in design weight. Effective 
skin thickness is based on a point design for the optimum skin/ 
stringer combination using a Lagrangian multiplier method (LMM) . 

A similar design approach is used to generate the optimum 0/+45/90 
composite laminates. This is a Level II program, though it does 
not look long. Only primary structural components are computed. 
Secondary weights are computed semi-empirically or empirically. 

The report contains information on the variation of material pro- 
perties with temperature and fatigue, which has been summarized 
in explicit form. The variety of construction types is representive 
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but limited to the available LMM derivations. A flutter design 
is included only as a lower bound to the stiffness via an opti- 
mality-criteria approach , though it considers only the first bend- 
ing and torsional modes. The point design procedure used is rel- 
atively sophisticated but does not address the combined synthesis 
problem. The program could possibly be used in a stand-alone mode 
with ACSYNT but would best be used as a data base for a number of 
technology items. The required input data are unknown. No cor- 
relations with existing aircraft were given. 

(2) The SWEEP program developed by North American Rockwell 
in reference 2 estimates complete vehicle weights. The fuselage 
and wing modules of interest here require up to 10 overlays and 
109 subroutines. The wing and empennage weight section has pre- 
liminary design to satisfy strength and stiffness. A separate 
module is available for specifying a flutter design obtained via 
optimality criteria as the lower bound for the strength design. 

The wing program is generally restricted to high-aspect-ratic 
wings with torque-box structures. Loads are generated from a 
station analysis though allowance for flexibility iterations is 
made. The design procedures employ direct numerical search or 
interpolation to find a strength design or to obtain the best 
stringer or rib spacings. The program contains a comprehensive 
data base for several types of construction (i.e., stringer, cor- 
rugated sheet, composites, sandwich, etc.). The weights of sec- 
ondary structural components such as flaps and ailerons are esti- 
mated from empirical data. 

The fuselage module contains detailed estimations for inter- 
nal and external geometry and loads for circular and a number of 
non-circular cross-sectional shapes applicable to both civil and 
military aircraft. Multilevel weight calculations are made for 
shell covers, major and minor frames, longerons, and bulkheads. 

In both modules weight correction factors for non-optimum weight 
are applied for a number of secondary weights. The program uses 
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multilevel types of analyses (input, rule-of-thumb, and detailed 
calculation) for secondary component weight estimation. Advanced 
structural technology includes some composite structures in the 
wing design. The use of optimality criteria in a companion flut- 
ter design program is the only apparent advanced numerical design 
technique. The portions of the code immediately usable by ACSYNT 
are the rule-of-thumb estimates of structural components and data 
associated with different construction types. Over two thousand 
inputs would have to be provided by ACSYNT to run the program as 
a direct module. Fuselage design takes between 10 and 60 seconds 
on a CDC 6600 computer. 

(3) The TSO program outlined in reference 23 is used for 
the aeroelastic design of isotropic and multilayered composite 
wing skin panels. The structure is modeled as a trapezoidal equiv 
alent flat plate, and stiffness and mass matrices are obtained by 
a Rayleigh-Ritz procedure. Optimization is carried out by sequen- 
tial unconstrained minimization techniques (SUMT) with a penalty 
function. Material distributions are described by continuous 
functional distributions. The program uses Woodward-Carmichael 
steady aerodynamics and kernel- function unsteady aerodynamics. 

No buckling criteria or transverse shear (spars or ribs) are in- 
cluded in the design. This is a Level II program. No secondary 
weight estimates are made. In general the program provides less 
detailed weight information than previously discussed. 

(4) The WIDOWAC (Wing Design Optimization with Aeroelastic 
Constraints) program outlined in references 24, 25, and 26 is 
currently being maintained at Langley Research Center for research 
studies in aeroelastic design. This is a Level II finite element 
structural program with symmetric shear web and membrane plate 
elements. It uses piston-theory supersonic aerodynamics and ker- 
nel-function subsonic aerodynamics. Optimization is by a SUMT 
search technique which uses approximate second derivatives in 
conjunction with Newton's method. Design-variable linking is used 
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to reduce the number of design variables. Run times typically 
require 1-5 minutes on the CDC 6600. Maximum stress, minimum gage, 
and flutter constraints are included in the design. No buckling 
or advanced material technology is included. No secondary struc- 
ture weight estimation is included. WIDOWAC should best be used 
in a stand-alone mode with ACSYNT for comparative purposes where 
flutter is of concern since no compatible automatic data genera- 
tion exists. 

(5) The SWIFT program developed at Langley is outlined in 
reference 27. This is a Level II program for the minimum-weight 
design of wings for combined strength and flutter requirements. 

The program uses an equivalent plate structural analysis and pis- 
ton-theory unsteady aerodynamics. No secondary structural weight 
is computed. Primary structural weight estimation is restricted 
to isotropic cover sheet material for a trapezoidal wing. No 
advanced structural concepts are used. This was a forerunner of 
the WADES program and has been superseded by enhanced versions. 

(6) The SSAM program was adapted for use in the ODIN synthe- 
sis program and is outlined in reference 28. The program is ap- 
plicable to the strength design of high-aspect-ratio swept wings. 

It uses the aeroelastic subsonic lifting-line theory developed by 
Gray and Schenk in reference 29. The box-beam structure is sized 
iteratively using a stress-ratio algorithm. Strength, buckling 
and minimum-gage constraints are included in the box structural 
model. No advanced materials such as composites are currently 
included. It may be used now as a stand-alone program with ACSYNT. 
Several non-optimum factors are included for the estimation of 
secondary structural weights. Automatic data generation would be 
required for use with ACSYNT directly. Execution times would be 

on the order of 1-2 seconds on a CDC 7600. The external loads are 
limited to the subsonic regime. 
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(7) The SAD program under development at Ames Research Cen- 
ter is a small in termed! ate- scale finite-element structural design 
program to support aerospace vehicle synthesis. This is a very 
efficient general structural design code for research in optimi- 
zation. It contains stress, displacement, and frequency con- 
straints. It has no plate buckling criteria, automatic data 
generation, or aerodynamic loads. Only primary structural weights 
are computed. No provisions for advanced structural concepts are 
currently incorporated. Superior advanced design concepts in- 
cluded are stress-ratio and feasible-directions numerical searches, 
inverse design space, Taylor-series constraint approximation, an- 
alytic gradients, design-variable linking, and multiple loading 
conditions. The only mode of operation with ACSYNT in the near 
future would be for stand-alone comparisons for variations in 
parameters. Configurations considered are limited only be genera- 
tion of geometric and load data. 

(8) The ACCESS program (refs. 30, 31, and 32) is a pilot 
program to study and demonstrate approximation-concept capabili- 
ties in the synthesis of general structures by means of the 
finite-element method of structural analysis. This is a very 
efficient Level II general purpose structural design code. It 
contains stress and displacement constraints). No plate buckling 
constraints, automatic data generation or aerodynamic loads are 
incorporated at this time. Only primary structural weight is 
computed. Flutter and composite materials are being incorporated. 
Advanced design concepts include NEWSUMT and feasible-directions 
numerical search techniques, design-variable linking, regionaliza- 
tion with respect to element configuration, inverse design varia- 
bles, constraint deletion, Taylor-series constraint approximation, 
analytic gradient information, and multiple loading conditions. 

The only direct usefulness in vehicle synthesis would be the 
generation of trend information in a stand-alone mode. The design 
concepts may be of use in various program developments. 
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(9) The ORACLE program used by Boeing for the preliminary 
design of high-aspect-ratio wings is documented in reference 33. 
This program uses the finite-element flexibility model outlined 
in reference 29 and also used by the SSAM program. This version 
uses more structural detail than used by SSAM in the beam model 
to size the structural box. Non-opt inrun weight fractions and 
secondary weights are also computed from statistical methods. 

Approximately the same data generation for use with ACSYNT would 
be required for this program as for SSAM. Demonstrated correla- 
tions exist in the above reference. 

(10) The APAS program was developed under Convair IRAD and 
is outlined in reference 34. The program can currently resize 
wing and fuselage components for strength requirements using a 
combination of beam and finite-element procedures. This is a 
Level II program primarily for transport aircraft. No documenta- 
tion of comparisons with existing aircraft was accessible for the 
prediction of non-optimum structural weight or for secondary struc- 
tural weights. A structural synthesis capability is available tor 
a wide variety of structural concepts including composite sandwich 
construction. A mixed optimization method such as reported by 
Sobieszczanski and Loendorf (ref. 35) was the basis of the design 
philosophy. The Fiacco-McCormick method using the Fletcher-Powell- 
Davidon unconstrained minimization technique was used as the design 
algorithm. The effects of multiple design conditions and fatigue 
can be considered. The geometry and external loads are input 
rather than calculated internally. The program appears to be too 
large for direct incorporation in vehicle synthesis. Documenta- 
tion would need further expansion to be of use as a source of 
information. 

A comparison of various program features and capabilities is 
given in Table II. This includes an assessment of the applicabil- 
ity of the program, the advanced technology included in the design, 
and the advanced design techniques used in sizing the structure. 
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The table haa been compiled from previously mentioned references 
and the information is dated. The programs examined to not rep- 
resent the complete list of programs available in industry. These 
programs were examined in depth primarily because of their acess- 
ibility or the availability of documentation. In general, large- 
scale finite-element programs were excluded because of their 
computational expense in a vehicle synthesis environment. 


TABLE II.- 

SUMMARY OF ANALYTICAL WEIGHT ESTIMATION PROGRAM CAPABILITIES 


B B P FE P B FE FE I B I B P 




W 

a 

Q 

< 

1 

i 


Program 

Capability 


ACSYNT Level 


Structural Model 1 

Applicability 
wings - High AR 
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Arbitrary Shape 

Automatic Geometry Definition 

Component Weight Estimation 
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Structural Technology 
Isotropic Materials 
Orthotropic Materials 
Composites 
Thermal Analysis 
Fatigue 


1 B * beam model, P * plate model, FE * finite element representation 
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TABLE II.- Continued 


SUMMARY OF ANALYTICAL WEIGHT ESTIMATION PROGRAM CAPABILITIES 


Program 


Capability 


ACSYNT Level 
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Loads 

Constant Pressure 
Piston Theory 
Modified Strip Loading 
F. E. Static Aerodynamics 
F. E. Unsteady Aerodynamics 
Discrete Loads 
Discrete Masses 

Aeroelasticity 

Flexible-Wing Loads 
Flutter 
Divergence 
Reversal 

Contro.. Feedback 

Structural Design Constraints 
Strength - Isotropic 

- Orthotropic 

- Composite 
Buckling - Isotripic 

- Orthotropic 

- Composite 
Interaction Curve 
Displacement 
Stiffness 
Frequency 

Flutter 
Divergence 
Minimum Gage 


Design Algorithm/Technique 
Feasible Directions 
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TABLE II.- Concluded 

SUMMARY OF ANALYTICAL WEIGHT ESTIMATION PROGRAM CAPABILITIES 


Program 

Capability 

FAIRCHILD 

SHEEP j 

TSO 

WIDOWAC 

SWIFT 

§ 

w 

SAD 

ACCESS 

3 

o 

APAS 

WADES 

ACSYNT Level 

II 

II 

II 

II 

II 

SI 

II 

II 

II 

II 

II 

Stress Ratio 


H 

■ 
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fl 

D 

■ 
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|| 

1 

i 

1: 
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1 
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II 
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■ 

| 
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II 
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■ 

Inverse Design Space 
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, Hr 


9 

If 

H 

D 

| 
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Lagrange Multiplier 


I 


SI 
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Optimality Criteria 

H 

1 

I 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

II 


A number of other programs available in the literature were 
examined. Their references and a brief description of salient 
features are included here. Most were not given detailed consid- 
eration because of their large-scale or proprietary nature. The 
othei programs examined were: 

ASOP - This is a general large-scale finite-element program 
using a modified stress-ratio approach for strength and a numeri- 
cal search for displacement constraints (refs. 36 and 37). A 
companion flutter and strength optimization program is described 
in reference 38. 

SAVES - This program uses NASTRAN, Carmichael-Woodward aero- 
dynamics, and stress resizing in an aeroelastic synthesis proce- 
dure for wings (ref. 39). 

FADES - This is a fuselage design code using finite-element 
analysis and mixed optimization in which partitioning into sub- 
structures is performed for resizing (ref. 35). 
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OPTIM II - This is a general-purpose finite-element program 
for minimum weight structures subjected to static loading condi- 
tions. Optimization is based on optimality criteria. It contains 
no aerodynamic interface (refs. 40 and 41). 

OPTSTATIC - This is a finite-element program for studying 
resizing statically loaded structures based on strain-energy- 
distribution optimality criteria (refs. 42 and 43) . 

OPTCOMPOSITE - This is a variation of the OPTSTATIC program 
adapted for the minimum-weight design of multilayered composites, 
(ref. 44) . 

NASTRAN - NASA-supported finite-clement structural analysis 
program (ref. 45) . 

SNAP - This is a large-scale finite-element program. 

ATLAS - This is a Boeing-produced large scale finite-element 
program for aircraft structures. It is a major aeroelastic anal- 
ysis program with a strength resizing capability. 

ASDP - This finite-element program is suitable for designing 
minimum-weight structures under static loading conditions. It 
uses a feasible-direction search technique (refs. 46 and 47). 

WINGOPT - This is the original program from which the WADES 
program evolved. It also includes geometric parameters of the 
wing for designing for optimum configurations (ref. 48) . 

ECI-ICES-STRUDL/DYNAL - This a large-scale finite-element 
structural analysis program. It uses design-table look-up for 
minimum-weight strength sizing (ref. 49) . 

ARROW - This is a McDonnell-Douglas program for automated 
design of large aerospace structures subject to static loading. 
Optimization by both nonlinear programming and optimality criteria 
is being used together with the large finite-element code, FORMAT 
(ref. 50). 
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Survey Summary 


The three classes of structural weight estimating programs 
examined were: I - Empirical weight -estimating procedures, II - 

Semi-empirically derived procedures, and III - Analytical design 
procedures. The detail and relative sophistication of the tech- 
niques increase in the order given. 

The Class I weight estimating programs provide the most rapid 
estimating technique. They are often as accurate within the 
bounds of their data as the more complex methods. This is espe- 
cially true when the estimation of secondary structure is in- 
volved. Only the WTSIZ program examined here provides subcompon- 
ent weight estimates. 

The Class II weight-estimating equations are generally based 
on obtaining an equivalent volume of material to satisfy a given 
load distribution. This allows for a moderate amount of flex- 
bility to compare technology used in sizing and to evaluate the 
sensitivity to major geometric variations but provides little 
consideration of changes in the critical loading or of effects 
of minor structural innovations. The statistical correlations 
required to compute secondary weights make these methods no more 
accurate than Class I equations. 

The Class III techniques offer the full scope of analysis 
and design capabilities. These methods have still depended on the 
use of Class I estimates of secondary weights to be of use in 
vehicle synthesis. Of the programs examined, a wide variation in 
methods and technology exists. The industry-developed programs 
(1,2,6,9,10)* contain better technology and estimation of secon- 
dary weight. 

These programs put considerable emphasis on detailed sizing 
of such structural components as stiffeners, ribs, rings, and 


Numbers refer to analytically based programs discussed in section 

1 . 3 . 
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bulkheads, and on the inclusions of all applicable constraints. 
This is consistent with the industry approach that the technology 
is more important than numerical methodology used to automate the 
desiqn procedure. The research-oriented programs (Nos. 3, 4, 5, 7, 

8) have demonstrated the applicability and superiority of numer- 
ical design techniques but are currently capable of providing 
only trend information. In general, structural detail has been 
sacrificed in order to obtain computational efficiency. 

Recommendations of Survey 

The recommendations presented here are based in part on the 
evaluations of the programs described and in part on the intended 
use of the results in vehicle synthesis. The accuracy of the 
results is based on the requirements for first or second level 
analysis detail. The program requirements will be discussed in 
terms of the estimation of wing and fuselage weights. 

Because of the optimization techniques used in vehicle syn- 
thesis many vehicle configurations must be examined rapidly. Ex- 
plicit continuous functions such as those derived from statisti- 
cal methods of weight estimation are best suited to vehicle 
design. The procedure recommended here is to generate analytical- 
ly based explicit functions from second level type structural 
weight estimating programs. This procedure would entail compu- 
tation of vehicle component weights by systematic variation of 
the vehicle design parameters about a nominal configuration. An 
explicit function developed from techniques as regression analy- 
sis would be derived for vehicle synthesis. 

The general characteristics that must be possessed by struc- 
tural wing weight estimating codes are that they model the re- 
quirements of the spectrum of configurations to be considered. 

For the transport and fighter aircraft expected to be examined 
the structural design code should be able to assess the influence 
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of wing flexibility on static loads and flutter; it should be able 
to model the skin as a stiffened metallic or composite plate and 
consider the effects of spar and rib spacing; and it should be 
able to estimate, at least by empirical methods, the secondciry 
structural weights. The preferable technique for analyzing the 
flutter margin for design should incorporate a finite element 
stiffness approach. For most medium to high aspect ratio wings 
a beam finite element should suffice. 

Of the second level programs considered the SWEEP program is 
the closest to being able to meet all the general requirements. 

The structural design for flutter and wing flexibility are nor- 
mally performed separately from the strength sizing, which though 
not preferred, would be acceptable. The ORACLE program would be 
acceptable for subsonic strength design but is proprietary and 
lacks flutter and supersonic aerodynamic capabilities. Both pro- 
grams have been correlated against wing weight estimates. Both 
programs would require consistent scaling of inputs to get proper 
sensitivity data. The WADES program with the stiffened plate 
skin and secondary weight equations should provide trend informa- 
tion for medium to low aspect ratio wings. The programs WIDCWAC, 
SAD, or ACCESS are still far from providing accurate total wing 
weight data, but they might be adapted to provide sensitivity in- 
formation for unusual configurations. They are best for providing 
new design methodology that should be useful in the future. 

The best compromise for a Level II structural design program 
for the estimation of fuselage weights is a modification of the 
SWEEP program to accept the geometry descriptors of the ACSYNT 
program. This should include a station analysis and design using 
a beam model. This model should be compatible with finite-element 
methods for .cer growth. The basic geometry within the program 
would not be significantly affected; only generation of certain 
parameters in terms of ACSYNT descriptors would then be necessary. 
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The secondary weight calculations wtuld be used directly. To 
obtain vehicle sensitivities a consistent set of scaling relation- 
ships for the inputs will have to be generated. 

Modifications to Wing Design Code 

In the phase I work of this contract a number of limitations 
in the program modeling were identified which, when corrected, 
should improve the comparisons with actual aircraft. The recod- 
ing of the WADES program was undertaken in the second phase of 
this contract to incorporate the recommendations. The modifica- 
tions incorporated and described here include: 1) expansion of 

the geometric planform definition to non-trapezodial shapes, 2) 
incorporation of a stiffened panel structural model, 3) adaptation 
of a trim loads calculation reacting the wing-body interface loads 
at the fuselage junction, 4) estimation of secondary structural 
component weights, and 5) development of an approximate flutter 
constraint. The descriptions of the theoretical derivations and 
implementations of the modeling and design changes follows. 

Geometric Planform Definition 

As a result of the initial correlation studies, it was 
recommended to modify the definition of the aerodynamic and struc- 
tural planforms to allow for a small number of planform discon- 
tinuities. This was considered necessary because in the compari- 
sons with actual aircraft the estimation of the weight showed a 
strong correlation with the planform area used to define it. The 
use of a single trapesoidal wing segment was therefore inadequate 
for a reasonable weight estimate. 

The original and the new planform description of the WADES 
wing model are shown in figures 1 and 2. In each case, both geo- 
metric descriptions are still acceptable inputs to the program. 
Figure 1 shows the two acceptable descriptions of the aerodynamic 
planform of the wing. The original description consisted of a 
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single trapezoidal region defined by the root chord, R, the semi- 
span, SPAN, and the leading- and trailing-edge angles, 6 t and 
0 2 . The new description consists of up to two adjacent trapezoid- 
al regions. This allows for the presence of one major planform 
break in the wing. The presence of an overlapping fuselage of 
width from the centerline, YFUS, is also included. This also will 
be used in the definition of loads to specify the point at which 
the wing loads are transferred into fuselage. Only the portion 
of the wing outboard of the fuselage will be considered to carry 
aerodynamic loads. This alternate description is defined by the 
leading-edge coordinates XSA and YSA, the chord length, CSA, and 
t/c ratio, TO, at each of NSA stations (NSA £ 3) . 

Figure 2 shows the two acceptable descriptions of the struc- 
tural planform of the wing. The original planform description 
consisted of a single trapezoidal region whose leading and trail- 
ing edges were defined as a fixed percentage of the local chord. 
The new description consists of up to three adjacent trapezoidal 
regions. The breaks in the structural planform defined here do 
not have to correspond to the aerodynamic planform breaks. The 
spanwise location at each of i stations (i <_ NSS) is defined 
as a function of its fraction of span by ETASS(i). Three values 
define the chordwise locations of the structural regions at each 
station. The first value, XCSS(l,i), is the x/c location of 
the leading edge of the structural region at the station. The 
second, XCSS(2,i), is the x/c location of the trailing edge 
inboard of that station. The third, XCSS(3,i), is the x/c lo- 
cation of the trailing edge of the structure outboard of that 
station. This arrangement requires that the leading edge of the 
structure be continuous. However, the trailing-edge structure 
may be discontinuous. At the root and at the tip, XCSS(2,i) 
equals SCSS(3,i). 

computation of wing weight components was modified to 
reflect the changes in planform definitions. Primary structural 
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weights will be computed only for structural regions defined in 
figure 2. The calculation of secondary structural components 
(i.e., flaps, ailerons, and leading-edge slats) has been incor- 
porated using empirical techniques. The details of the empirical 
equations that have been incorporated are given later in this 
report. 


STIFFENED PANEL STRUCTURAL MODEL 

Modifications to the WADES program during this period were 
directed towards implementing the proposed changes in the struc- 
tural model to include the effects of stiffening and buckling of 
the wing covers. In order to implement the modeling of stiffen- 
ers an approximate "Z" stiffened plate description was developed. 
The variables used to describe this plate are discussed below. 

The original description of the cover sheets of conventional 
wing structures used by the WADES program consisted of a single 
distributed function, T(5,n). It was assumed that this was 
equivalent to the stiffened plate in total volume of material at 
the optimum design. It did not include its orthotropic effects 
on the wing stiffness. The initial description chosen to repre- 
sent the skin-stiffener arrangement is shown in figure 3. The 
wing segment in figure 3 is partitioned into chordwise and span- 
wise panels. For the purposes of simplicity the ribs are assumed 
to be parallel to the x-axis , with a spacing in the spanwise 
direction of YRIB. This is typical of medium to low aspect ra- 
tios and sweeps. The structural reference angle, THET, has been 
included to define the principal axis along which the structure 
will be analyzed. This axis should be parallel to the primary 
direction of the load due to bending. Though some allowance for 
taper in the width of the stiffened panel is permitted, the in- 
ternal loads will be estimated only parallel to the structural 
axes . 
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The panel cross-section in figure 3 shows the variables used 
to describe the skin and stringer dimensions. An integral "Z" 
stiffener is shown in this derivation, though the model is not 
restricted to that configuration. The isotropic skin thickness is 
again represented by T(£,n). The variable describing the stiff- 
ener spacing is Bl. The stiffener dimensions are the width of 
the web, B2, and the flange width, B3. For simplicity the flange 
and web thicknesses are both equal to T2 . The thickness of the 
spar is shown to be T3 with a semi -height of D. 

In order to facilitate the use of the variables in the de- 
sign procedure each of these variables has been approximated by 
a function. This is necessary in order to incorporate them in 
the continuous function analysis approach used by the WADES pro- 
gram. The design variables are then the coefficients of the 
appropriate functions. This linking in design variables also 
reduces the total number of variables required. The stiffened 
plate is thus represented by the thickness function, T(£,n), the 
stringer thickness, T2(n), and stringer dimensions, B2(n) and 
B3(n). The stringer spacing, however, is represented by the in- 
verse of the spacing, as Bll(n) = 1/Bl . This design variable 
was chosen in order to allow explicit integration in the weight 
equations. The stringer dimensions have initially been repre- 
sented as functions of n only. This was an arbitrary simplifi- 
cation, and the dimensions may be expanded to functions of C 
and n if required. The code is being arranged so that any of 
the stringer variables may optionally be left in or left out of 
the design. 


STRUCTURAL DESIGN 

The variables used to describe and analyze the stiffened 
plate were described above and in figure 3. Their incorporation 
in the analysis and design is now described. 
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Where the inplane stiffness of the cover plates was described 
by the product of the modulus and the skin thickness , T, in the 
old model, it is now obtained a6 the superposition of the stiff- 
ness of the skin and that of the stringers. The stiffness is then 
written as 

(E»TJ - (Eskinlt ♦ [R) T [E gtr ] [R] (B2 + 2 B3)T2 BIl (1) 


where (RJ is transformation rotating the modulus of the string- 
er along the structural axes. 


The weight of the skin and stringers may be obtained as the 
sum of the weight of the skin, w S ki n ' an£ * t * ie we i 9 ^t of the 
stringers, W . They may be obtained by direct integration cf 
the following equations: 


W 


skin 


p t (C.n) dA 
1 A 


( 2 ) 


and 

W . * p (B2 + 2 B3) T2 BIl dA (3) 

str Ja 

Here A is the structural planform over which the integration 
is carried out, and p is the material density. These weights 
are also used as the objective function during optimization. 

The design procedure used to size the wing for strength can 
no longer be linearized as before. The procedure used now will 
be similar to that used Lo size the composite design. Instead 
of the von Mises' stress being held constant during a sizing 
cycle, the component edge loading will be held stationary. This 
technique still allows for changes in load path between skin and 
stringers during the optimization sequence. Thus, the addition 
of the buckling imposes a new nonlinearity in the constrr ! ■ _s. 

The strength constraints on the wing design that will be 
considered are the von Mises' stress resultant in the skin, the 
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panel buckling of the skin panels, the buckling of the stringers, 
and the minimum gage requirements on both. The panel buckling 
will be in the form of a standard interaction formula for combined 
shear and compressive buckling of a rectangular panel. Local 
buckling of the panel will be calculated assuming the lower bound 
buckling coefficients for panels of infinite aspect ratio. The 
strength design must satisfy the following inequalities: 


Von Mises’ stress: 

0 lrM < S 
VM — max„ 
c 


(4) 


Buckling interaction: 

1 i R C * *5 


£ 

c 

F 

max 



(5) 


Here f c is the compressive stress along the structural axis and 
f is the shear stress. The buckling allowables are obtained as 
averages of the tensile and compressive allowables. This is used 
to average the difference in gages of the upper and lower skins. 
The longitudiral and shear allowables may then be written as 


F max * l ' 2 


“* 

f4TT 2 E„ (T/Bl) 2 


S + min 

. S 


max fc 

10/1 2 \ max^ 

[ 12(1 “ \i ) cj 



( 6 ) 


and 


F = min 
scr 


5.62TT 2 E r (T/Bl) 


12(1 - y 2 ) 


max 


(7) 


where E r is the reduced modulus of the material, and 

, and t are the ultimate compressive, tensile and shear 

max^ max F 

allowables, respectively,, 

Similarly, the stringers are sized by their allowable buck- 
ling stresses. No shear effects are included in the sizing of 


max. 
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the stringers 


The allowable stress for the web is written as 


web 


1/2 


S v + min 
max t 


S , 
max c 


4n 2 E r (T2/B2) 2 ) 


12(1 - u 2 ) 

and the allowable for the stringer flange is written as 


F flange ' 1/2 



r 0 . 426 it 2 E (T2/B3) 2 

" 

S „ + min 

S , E 


max fc 

max c 12 (1 - y 2 ) 

*■ 


Trim Load* Calculation 

In the previous version of WADES, the externally applied 
lo«. ds were calculated from one of three methods: piston theory, 

constant pressure wing loading, or modified strip loading. The 
appropriate method was selected at program load time and was 
limited to that method for the remainder of execution. The pre- 
sent program version was modified to incorporate two types of 
load calculations - modified strip loading and either piston 
theory or constant pressure wing loading. It is intended that 
these two methods would represent subsonic and supersonic aero- 
dynamic loads. The choice of the method is specified by the user 
through the analysis option control parameter, IANAL (2 , IFLT) . 

The choice of piston theory or constant pressure wing loading 
routines is controlled at program load time. Only one of the 
methods may be loaded during a given run. 

Longitudinal trim of the aircraft was incorporated into the 
loads estimation in order to assess the weight penalty due to 
center of gravity travel. The center of gravity location is spec 
ified as the fraction of the mean aerodynamic chord forward of 
aerodynamic center (PMAC) . The equivalent tail load (P T ) and 
gross lift (GLR) required to trim the aircraft are then computed. 
Figure 4 shows the relative location of the aircraft center of 
gravity at the stress gross weight (SGW) , the mean aerodynamic 
chord and the tail load. The gross lift required and the 
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tail load are computed from the following summation of forces and 
moments: 


GLR - N,g *SGW ♦ P T ♦ ZF i - 0 


GLR-X^c - V SGW0 W - PMAC C^) + P^ + - 0 


Then the gross lift required becomes: 


GLF 


V SGW(X MAC - PMAC<C MAC - X HT } - LF i X i j X HT LF j 

X MAC " X HT 


( 10 ) 


and the tail load required for trim is: 


P T * -GLR + N Z «SGW - EF i (11) 

The reaction of tail and fuselage loads into the wing is 
incorporated as an equivalent force couple at the wing-body junc- 
tion. That is, the total force and moment of the body and tail 
loads are resolved into two forces at the leading and trailing 
edges of the structural planform. Figure 5 depicts the relative 
locations of the two equivalent reactions between the wing and 
fuselage taken from figure 4 at the y-location of the wing-fuse- 
lage junction (¥p US ) . The reactions are then included in the 
structural analysis as concentrated forces. The reaction forces 
are computed from the summation for forces and moments as follows: 

~ N z’ W BODY (X B~ X R 2 ) + P T (X HT _ ^2 } v 

R 1 V U2 ’ 


*2 “ f'-V W BODY + 


P T 


R ll 


(13) 
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The calculated load distribution on the wing varies slightly 
between the modified strip loading and the constant pressure or 
piston theory loads. In the first, the input values of the lo- 
cation and length of the mean aerodynamic chord are used only as 
a first guess. The values are subsequently recalculated and the 
updated values used in later calculations. In the second, the 
trim conditions are dependent on the input values of 


Critical Load Profile 

The identification and specification of the critical load- 
ing conditions to be used to size the primary structural ele- 
ments of the structure form a necessary part of determining the 
design loads. Here it is intended only to provide guidelines tc 
assist in that selection of flight conditions. 

The basic strength requirements must satisfy PAR Part 2b 
for civil aircraft and MIL-A-008860A or its equivalent for mil- 
itary aircraft. In general the strength requirements must be 
met for every combination of velocity and load factor within the 
maneuvering and gust envelopes in figure 6. Typically, the max- 
imum level speed (V c ) and the design dive speed (V^) are the 
critical gust conditions and their determination is the primary 
concern. Figure 7 is a plot of a typical structural design air- 
speed profile versus altitude. Typically the dive speed or dive 
Mach limit (M^) can be specified in terms of V c or the equiva-’ 
lent cruise Mach limit (M c ) , so that only V c or M Q must be 
determined. 

Once the speed profile of the mission is determined the 
three remaining parameters to be specified are the maximum ma- 
neuvering load factor (n ) , the maximum gust load factor (n ) , 

z y 

and the structural gross weight (SGW) . In lieu of special re- 
quirements suggested values for the limit maneuver load factors 
from MIL-A-8861 are as follows: 
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limit maneuver load factor - n 2 


Vehicle 'Class 

subsonic 

supersonic 

negative 

flaps 

Fighter, Attack 

8.0 

6.5 

-3.0 

4.0 

Bomber T. 

4.0 

4.0 

1 

to 

• 

o 

2.5 

Bomber II 

3.0 

3.0 

1 

H* 

• 

O 

2.0 

Cargo Assault 

3.0 

3.0 

1 

• 

o 

2.0 

Cargo Transport 

2.5 

2.5 

1 

• 

O 

2.0 


The limit gust lead factor may be calculated either from a 
theoretical gust profile or from the appropriate FAR or MIL em- 
pirical equation. The FAR-25 equation for the estimation of 
the gust load factor coded in the WADES program is as follows: 


n z = 1 


K g U de V a 
498 (W/S) 


(14) 


where 


0 . 88p 

K_ = t — rr- ;° = gust alleviation factor 


g 5 . 3+p 


(15) 


p - 2 (W/S) _ a i r pi ane mass ratio; 

9 p cag 


(16) 


= derived gust velocities (fps) ; 

p = density of air (slugs/cu. ft . ) ; 

W/S = wing loading (psf ) ; 

C = mean geometric chord (ft.); 

g = acceleration due to gravity (ft./sec.); \ 

V = airplane equivalent speed (knots) ; 

a = slope of the airplane normal force 
coefficient curve Cn^ per radian 
if the loads are applied to the wings 
and horizontal tail surfaces simultan- 
eously in a rational method. f 



The derived gust velocity, U^, is also specified by FAR-25 
to be computed from the following: 

f rough-air gust at V 


[66 - 


66 fps 


h-20000 

30000 


2 8] fps 


gust at V 


50 - 


50 fps 


13-20000 1 f 

30000 1 J *■ 


gust at V 


D 


25 fps 


^ c . h-20000 t t C 1 

25 + ' 30000 12 - 5 J £ P S 

L 


0 < h < 20000 


20000 < h < 50000 


0 < h < 20000 


20000 < h < 50000 


0 < h < 20000 


20000 < h < 50000 (17) 




The structural gross weight used in the estimation of the 
critical maneuver loading condition is the weight which produces 
the maximum stresses in the wing or produces the maximum stres- 
ses in the wing or produces the maximum load on the fuselage. 

The cases which should be checked for this condition include: 


1) For fuel carried in the wing, the maximum vehicle weight 
with only the fuel reserves left in the wing. 

2) The maximum vehicle weight with zero fuel in the wings. 
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3) The maximum weight at which the maximum maneuver load 
factor may first be achieved as specified by the mission profile. 

The additional loading conditions, which may be critical in 
the flight and ground handling envelope, are as follows: 

1) The maximum load factor in the extended flaps conditions 
during takeoff. 

2) The maximum aileron roll condition with flaps retracted 
and extended during takeoff. 

3) The lg trimmed landing impact condition. 

4) The maximum braking roll condition. 

SECONDARY WEIGHT EQUATIONS 

In order to estimate the weight of leading- and trailing- 
edqe devices, ailerons and fixed secondary structure, a set of 
empirical equations was included in the WADES program. Refer- 
ences 2 and 22 were examined for possible equations. Because of 
the additional detail and documentation included, the equations 
of reference 2 were added to the WADES program. The use of these 
equations does i.aquire some additional user input. Appropriate 
values are calculated from the WADES geometry description where- 
ever possible. The calculated program values may optionally be 
over-written by the user where specific secondary structural 
unit weights are known. 

The cipproach taken in the SWEEP program (ref. 2) was to es- 
timate each major leading- and trailing-edge component with 
statistical equations based on component geometry parameters 
and/or vehicle design criteria. The basic equations for control 
surface devices are modified so that the unit weights can be 
adjusted through specific types of data in the input. The form 
of the equations has been derived using a general equation form. 

The weight estimation equation can be expressed in general form as: 
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( 18 ) 


where 


w/s * K iKj + IAK] ( (w/s) ) 


w/s 

K 

o 

K » 

IAK 

(w/s) o 


estimated unit weight 

general weight coefficient to be 
used by the user, 1.0 unless changed 

basic statistical equation correlation 
factor, different for all components 

derived unit weight modification factors 

basic statistical unit weight (function 
of vehicle environment and geometry) 


For all devices, provisions are made in the input data set and 
analysis logic to allow the user to specify desired unit weights, 
in lieu of the program derived data. 

The basic statistical unit weight, (w/s) Q , is derived for 
each component in one of the following three forms: 


w/s = C t (C 2 X, + C 3 ) 


( 19 ) 


w/s = C (X ,) Cj + C (X ) C7 

4 Z b 3 


( 20 ) 


w/s = C e (XJ 9 


( 21 ) 


where the Cj_ 9 a * e equation constants, and the are 

estimation parameters based on vehicle criteria and component 
geometry . 

The unit weight modification factor IAK consists of three 
terms : 


IAK = AK, 


(t/c) 


ref 


(t/c) , 


0.25 


r 0*125 1 

+ ak 3 + ak._ In •• i.oj 


( 22 ) 
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where 


AK = basic incremental factor for thickness ratio 

2 

AK * basic incremental factor for available volume 

3 

AK * basic incremental factor for number of actuators 

4 

(t/c) , = 0.10, constant 

' ref 

(t/c) . = aerodynamic thickness ratio at midspan of each 

1 segment panel 

N = number of actuators per segment panel 

Term 1 of equation (18) is included in the fixed structure equa- 
tions . 


Leading-Edge Structure 

The leading-edge structure is assumed to include all struc- 
tures forward of the front spar and between the reference lines 
defined by the y-coordinate of the wing-fuselage junction and 
the outboard tip of the structure. Four leading-edge unit weight 
equations for various types of devices have been included and 
are optionally available through the program parameter, ISEC(l). 
The optional unit weight equations ares 

ISEC(l) = 0, user inputs (w/s) . 

13EC ( 1 ) = 1, fixed leading-edge structure: 


(w/s) = 0.00077 


0.8 Q S. 
max le 


ave 


+ 0.83 


(23) 


( l - 


(w/s) w = K w \ 1.50 + 0.10 


0 . 2 5 

0.10 1 \ 

( t/c> ave j i 


(w/s) (24) 


where 

^max = max ^ mum dynamic pressure, generally determined at 
V , sea level 
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S le “ ex P ose< ^ theoretical leading-edge planform area 

C * average chord determined by dividing exposed area 
by exposed leading-edge span measured along front 
spar 

(t/c) * aerodynamic thickness ratio at exposed leading-edge 

ave midspan 


ISEC(l) = 2, leading-edge slats: 


(w/s) Q ■ 


0.551 


f >,0.32 

N ,.DGW 
zult 


w 


+ 1.0 


0.8Q S . 
max pnl 


pnl 


0 2 S' 


(2 


(w/s) sl = K sl< 1,0 + 0,10 


0.10 


(t/c) 


avej 


0 . 2 S 


+ 0.01 + 1.0 


0.125 

N 


- i.on 


(w/s) 


(2 


J 


ISEC(l) = 3, leading-edge Kruger flaps: 

0 . 3 2 


(w/s) o = 


0.413 


N 


zult DGW 


w 


+ 0.667 


0 . 8Q S . 
max pnl 


pnl 


0 2 5 


(2 


(w/s ) kr = K kr < !.° + 0.10 


0.10 


(t/c) 


ave 


0.25 


+ 0.01 + 0.75 


0 12 5 

N * - 1.0 


(w/s) 


(2 


I SEC ( 1 ) = 4, droop leading edge: 


(w/s ) q = 

0.00077 

r ° . 8Q msav s „„ , ] 

max pnl 

c 

+ 0.83 

+ 

0.33 



ave J 

J 




0 . 8Q S . 
max pnl 


pnl 


(2 
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- K dn( 1 - 725 * 0 - 10 LffiJ *'” * °-° l 


+ 0.50 N°* 125 - 1 


■•]} 


(w/s). 


! z ul t = u ^ t -i mate positive load factor 

)WG = basic flight design gross weight 

> w = gross wing planform area 

= planform area for each device segment 

> . = device segment span measured along forward 

pn device control line 

: = average device segment chord 


Trailing-Edge Structure 

Trailing-edge structure unit weights may be computed for 
three different types of devices: flaps, ailerons, and fixed 

trailing-edge structure. The effects of multisegment flap de- 
vices are included in the correlation coefficient, K t yp e * All 
trailing-edge structural planform not specified as either flap 
or aileron is considered to be fixed trailing edge. The calcula- 
tion of the appropriate device is controlled by the option para- 
meter, ISEC. The following unit weight equations have been 
incorporated . 

Trailing-edge flaps - ISEC(2) = 0, user inputs (w/s)* 

ISEC (2 ) > 0, 


(w/s) = 0.69 

o 


14 -< <W b f 

100 (t/c) 


0 • i 5 


(w/s) 


f * K f K type 


0.01 + 1.5 N°* 125 - lj^ (w/s) 


[i - ISEC (2 ) ] 


( 32 ) 


Fixed trailing-edge structure (secondary structure) - the basic 
statistical equation correlation factor, in the general 

equation 18, is adjusted for fixed trailing structures by a 
coefficient that is sensitive to the maximum design dynamic 


pressure. The correction factor AK^ is determined as: 


AK = C 
q a 



®max 

v; 





(33) 


Here C a and are constants, currently assigned values of 


1.0 and 0.70 for wing and 0.75 and 0.70 for horizontal and ver- 
tical tail surfaces, and 


Q o = reference dynamic pressure, 950 psf 


^max * max ^ muin dynamic pressure, psf 


ISEC(3) = 0, user inputs (w/s), 
ISEC (3) = 1, 


(w/s) Q = 0.0165 


0.35 Q 

max te 


te 


+ 1.45 


(34) 


(w/s) w = K w ^ 1.0 + AK g + 0.10 


0.10 


(t/c) 


ave, 


0.25 


(w/s) 


(35) 


Ailerons - ISEC (4) = 0, user inputs (w/s) 
ISEC ( 4 ) = 1, 
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..-J* 


(w/s) Q 


^0 .01825 - a j + 1.55 + 0.50 0* 35(Q max ) 



(w/s) a 


K a ^i • ° +0.10 

[ 1 ^ 

N°* 125 - ij \ (W/S) Q 


0.10 


(t/c) 


ave 


+ 0.01 
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The K type f actor in equation 32 is selected from a table 
of factors based on the type of flap specified. The following 
table values for K type a i° n 9 with the indicator control word 
values . 


FLAP-TYPE INDICATOR AND CORRELATION COEFFICIENTS 


Flap Type 

Indicator 
Value 
ISEC (2 ) 

Correlation 

Coefficient 

K 

type 

Simple 

1 

1.000 

Single-slotted 

2 

1.250 

Double-s lotted 

3 

1.500 

Triple-slotted 

4 

1.750 


The flap segment area, S^, found in equation (31), is the 
sum of all chordwise panel areas. Thus, for triple-slotted 
flaps, actual planform areas are computed for each of the three 
chordwise panels. is then the sum. 

FLUTTER SENSITIVITY 

In the original version of the WADES program that was inte 
grated with ACSYNT the design loop that performed simultaneous 

4 : 





strength and flutter optimization was left out due to core limi- 
tations. In order to get flutter back into the synthesis loop, 
it is anticipated that an approximation concept flutter constraint 
can be derived. In order to test its feasibility for use in the 
WADES program a sensitivity run was attempted to compare a flut- 
ter approximation with actual flutter computations. In this 
example an initial strength design was first obtained and then 
each design variable was perturbed about that design point. The 
flutter speed was evaluated for each perturbed design variable. 

The computed flutter point was then compared with a Taylor series 
approximation of the flutter dynamic pressure obtained by finite 
difference. Figure 8 is a plot of two such sensitivity runs for 
two of the ten design variables used in the example. The flutter 
dynamic pressure and its Taylor series approximation are plotted 
versus wing weight. Because the design variables are coefficients 
of a polynomial function the perturbed variables were selected 
so as to obtain a unit change in the wing weight. The range of 
given variable changes spans about a 25% change in the total 
wing weight. 


Flutter Constraint Approximation 

A constraint based on the sensitivity information discussed 
in section 2.8 was incorporated in the WADES program. The basic 
procedure used was: (1) analyze initial structure; (2) generate 

the gradients of the flutter speed by finite difference; (3) ex- 
press the flutter dynamic pressure, q^, and Mach number, Mf, as 
Taylor series expansions in terms of the design variables, x^; 
and (4) include the series in the following constraints on the 
flutter dynamic pressure and Mach number in the strength design 
iteration: 

Minimum flutter dynamic pressure: 
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Minimum flutter dynamic pressure: 


NDV 3M, 


M 


fo, 




= 1 


E 

i=l 


■SxJ- (x i ‘ V 


FMMIN 


< 0 


(39) 


k=l, NOFC 


where 

QFMIN = minimum allowable flutter dynamic pressure 

FMMIN = minimum allowable flutter Mach number 
( ) = initial value of parameter 

After the structure is resized, the wing weight is checked 
for convergence and the design procedure is either restarted or 
terminated if no change was observed during the last iteration. 
The gradient calculation may be computed every iteration or may 
be updated periodically in order to reduce computational effort. 
The origin, about which the Taylo*' series is computed, is up- 
dated during ev* ry iteration. 

A preliminary example of this procedure for flutter design 
was executed using the combined strength and flutter wing design 
problem cf Stroud. Ten cycles of the strength design were per- 
formed to obtain an initial material distribution that satisfied 
the strength constraints but violated the flutter Mach number 
requirement by 10 percent. Two iterations of combined flutter 
and strength optimization were performed. The reanalysis of the 
designs using an approximate flutter constraint after each 




iteration showed that the actual flutter Mach numbers were within 
0.4 percent of the estimated values. 

The user sheets and descriptions were also updated to re- 
flect the new program variables. The program descriptions were 
expanded to include the new planform definitions/ buckling, and 
secondary weight calculations. A few new variables will be re- 
quired to execute this program version and will be noted in the 
documentation. The changes in old variable descriptions or 
definitions will also be noted. 

Comparison of WADES Estimates with the F-5A Wing 

In phase I of this contract a number of areas for improve- 
ment were identified and changes recommended. The previous sec- 
tion described their integration into the WADES program. The 
computer code was used to recompute the estimates of the F-5A 
wing weights performed in phase I. The comparisons of the ma- 
terial distributions, estimated component weights, and spanwise 
loading distribution computed by the WADES program follow. The 
results shown here reflect the modifications described in the 
previous section. A comparison with previous estimates is not 
given here but may be obtained from reference 1. 

Structural Model of F-5A 

A re-correlation of the F-5A weight calculations was per- 
formed. The input data were modified to reflect the changes in 
required inputs. The changes were made primarily to reflect the 
program changes in planform description of the aerodynamic and 
structural models, in the calculations of loads for trimmed 
flight, and in the addition of the calculation of secondary 
structural component weights. The use of the multiregion plan- 
form descriptions of the various portions of the wing is described. 

Figure 9 is a pictorial representation of the geometry of 
the mathematical model used by the WADES program to analyze the 
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F-5A wing. The regions outlined in the interior of the wing 
planform show the use of the new structural planform description. 
The three segment structural planform follows the carry-thru, 
intermediate and outboard panels used in the F-5A as primary load 
carrying structure. The leading edge and trailing edge flaps 
and the ailerons are indicated where secondary structural com- 
ponent weights have been computed for these devices. The actual 
dimensions of these devices have been used in the empirical 
weight calculations. The remaining planform is designated as 
secondary structure, and a unit weight estimated as a fixed 
trailing edge structure. The planform area used in the estima- 
tion of aerodynamic loads is that area outboard of the wing-body 
intersection. Structurally, only that planform area is consid- 
ered to carry aerodynamic loads. No weights are computed for 
the planform area interior to the vehicle fuselage except for 
the wing carry-through structure. 

Within the structural planform in figure 9 is a contour plot 
of the functional representation of the F-5A skin thickness dis- 
tribution generated from the upper surface of the wing cover. 

This is the functional representation, t(£,n)» used by the WADES 
program to analyze the F-5A. A least-squares functional fit of 
the material distribution over the surface of the wing in refer- 
ence 51 was computed for both the skin thickness, t(£,n)# and 
the equivalent thickness, t(£,n), in order to compare the analysis 
results predicted by the WADES program. The equivalent thickness 
was computed by distributing the additional spar cap material as 
wing covering in order to assess the inclusion of the bending 
effectiveness of the spars. These functional approximations of 
the skin thickness are used to check the accuracy of the analysis 
model and as initial values to the design procedure. Subsequent 
designs are compared with these to assess their validity. 
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Distribution of Material 

Four cases were used to check the estimation of the struc- 
tural material distributions by the WADES program. The first 
two analyze the F-5A wing using the thickness functions obtained 
from the surface fits of the skin covering, t(£,n), and the 
equivalent thickness, t(£,n). Only the static analysis was per- 
formed for these cases. The third and fourth cases used the 
first two material distributions as initial functions and rede- 
signed the material distributions to satisfy the constraints 
defined by the program for sizing only the skin thickness, t(£,n)» 
and for sizing the skin thickness with equivalent spar caps. 

The spar caps were approximated by including a stringer of mini- 
mum gage. The material gates were resized to satisfy the buck- 
ling, strength, and minimum gage requirements for the stresses 
generated by the WADES program. The estimated of these four 
examples were then compared to data obtained from the F-5A. 

Figure 10 is a plot of the cross-sectional area of struc- 
tural material versus span for the F-5A and the four WADES check 
cases. The F-5A structural areas were those used in reference 
52 to determine the margins of safety for the wing stress analy- 
sis. The cross-sectional areas estimated by the WADES program 
were obtained by integrating the skin thickness parallel to the 
structural reference axis along the 35 percent chord. The dis- 
continuities in the material distributions are due to discrete 
changes in the planform describing the structure. The estimates 
by the WADES program are generally within 15 percent of the actual 
distributions. In the first two cases, the functional distribu- 
tions obtained from surface fits are below F-5A net areas over 
most of the span. The low values indicate that a certain amount 
of additional material that contributes to the bending strength 
is not included in the equation fits. The structural areas es- 
timated by the sizing a.’ orithm are mixed both high and low in 
various segments of the wing. Some material is missed in the 


skin resizing cycle at wing stations outboard of WS 100 due to 
insufficient loading conditions. This region is sized by mini- 
mum gage and external store separation criteria. 1 ' low area 
value at the root for the redesign of t(5,n) is part due to 
clamped boundary condition and due to the over-estimation of the 
depth of the carry-through structure by the depth function. 

Figure 11 is a plot of the estimated moment of inertia , lx, 
of structural material versus span for the F-5A and the four 
WADES check cases. The F-5A structural areas were eff active 
moments of inertia used in the wing stress analysis in reference 
52. The moments of inertia estimated by the WADES program were 
obtained by integrating the product of the skin thickness and 
the square of the wing depth parallel to the structural reference 
axis. The inertias generally follow the trends displayed by the 
areas. The excessive inertias through the carry-through struc- 
ture are primarily the result of an over-estimation of the depth 
of the structure at the root. The depth function was derived to 
approximate the theoretical t/c at the root, rather than the 
constant section carry-through structure. 

Table III compares the weights predicted by the WADES pro- 
gram for each of the four check cases to the values obtained 
from the F-5A wing group weight statement. Of the six component 
items in the table estimated by the WADES program the center and 
outboard section weights are computed from structural analysis. 
The secondary structural weights, ailerons, and flaps are esti- 
mated from empirical techniques. For the purposes of generating 
weight correlation factors the center and outboard sections and 
the secondary structure are combined as one parameter for the 
basic structural weight. At the bottom of the table are the 
ratios of actual to computed groups. The first is the ratio of 
the net weights. The second is the ratio of the sum of the 
center, outboard, and secondary structure weights. The third is 
the ratio of the sum of the aileron and leading and trailing 
edge flap weights. 
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Of the sample cases in the table , the functional fit cf the 
skin thickness , t, underestimated the net weight , the available 
structural area, and tha moment cf inertia. The functional fit 
of the effective skin thickness, t, and the thickness distribu- 
tion from the redesign of t gave the most consistent comparison 
with F-5A areas and inertias and their predicted weights were 
within one percent. The use of the stringers in the redesign 
procedure to represent spar caps predicted the best weight, but 
over-estimated the structural areas and inertias. The average 
ratio of actual to computed total weight of the four check cases 
was 1.23, As a group the aileron and flap weights were within 
nine percent of the actual values. 


TABLE III 

Comparison of Estimated and Actual 
Group Weights for the F-5A 





Skin function fit 

Redesian 


Wing Component 

F-5A 

t 

t 

t 

r 

1 . 

Center Section-Basic Structure 

122.3 

159.7 

170.6 

157.3 

182.3 

2, 

Outer Section-Basic Structure 

716.6 

360.8 

411.1 

415.3 

446.2 

3. 

Secondary Structure 

36.6 

119.6 

119.6 

119.6 

119.6 

4. 

Ailerons 

35.7 

38.3 

38.3 

38.3 

38.3 

5. 

Flaps - Trailing Edge 

61.2 

51.4 

51.4 

51.4 

51.4 

6. 

- Leading Edge 

69.9 

63.6 

63.6 

63.6 

63.6 

Total 

1041.7 

793.4 

854.5 

845.7 

901.4 


(W ,/W ) 

act comp 

1.00 

1.31 

1.22 

1.23 

1.16 


(W /W ) , „ „ 

act comp 1+2+3 

1.00 

1.37 

1.25 

1.27 1 

1.17 


(W VW ) . c r 

act comp 4+5+6 

1.00 

1.09 

1.09 

1.09 ' 

1.09 


Externally Applied Loads 

Two modifications previously described have significant 
impact on the agreement of the estimated externally applied loads 
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with F-5A shear and moment distributions. The inclusion of the 
tail load required for trimmed flight increases the total lift 
on the wing, and the reaction of the resulting body weight at 
the wing-fuselage junction pushes the shear and moment distribu- 
tions outboard. Figures 12 and 13 are plots of the limit span- 
wise wing loads versus span for the F-5A and the loads estimated 
by the WADES program. Figure 12 shows the shear and moment dis- 
tributions for flight condition 123C-5, one of the maximum sym- 
metric loadings which designs the F-5A. Figure 13 shows the 
loading distributions for the dynamic landing condition, 

358B t - 118. 

The F-5A was structurally sized by three flight conditions 
in these check cases. The maximum symmetric load factor flight 
conditions, 104 and 123C-5, and the dynamic landing case, 358, 
were modeled as the critical loading cases in WADES. The re- 
sulting analysis and design indicated that the maximum symmetric 
load factor conditions sized the WADES design. The loads pre- 
dicted by the WADES program for condition 123C-5, in figure 12 
were within 5-10 percent of actual values for the shear, Vz , 
and bending moment, Mx, over the entire span. The torsional 
moment, My, was about 20 percent low at the root. However, no 
pitching moment was input into the aerodynamic computations . 

The dynamic landing case shows fair correlation at the root but 
varies considerably over the span. The Northrop loads were 
generated from the dynamics of the landing profile; the WADES 
loads were estimated as a static equivalent load at impact simu- 
lating the force required to absorb the kinetic energy due to 
the aircraft sink rate. 

Estimation of Fuselage Weights 

Development of an analytically based computer code for the 
estimation of aircraft fuselage weights was identified as the 
second contract item. The derivation of this code was originally 
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intended to be developed with the geometry descriptors of the 
ACSYNT program and to be called directly as a module of ACSYNT. 

The development of a fuselage design code was intended to pro- 
vide a means of assessing analytically the effects of advanced 
technology in aircraft design. Part of the survey of available 
computer codes given earlier in this report was directed towards 
assessing aircraft fuselage design programs to satisfy these 
requirements . 

Three factors were considered in establishing the approach 
taken in the development of the fuselage code: the recommenda- 

tions of the computer code survey, the requirements of the ACSYNT 
program, and the development of other codes for structural design 
to complement the ACSYNT program. 

The original recommendations of the survey of available 
computer codes were directed towards assessing the programs' 
potential for rapid estimates of structural weight or Level I 
weights estimation. In these terms the survey suggests that a 
combination of the procedures outlined by Shanley (ref. 15) and 
methods employed by the SWEEP program (ref. 2) should be adapted. 
The basic structural components would be estimated analytically 
with a great number of secondary structural component weights 
adapted from the SWEEP program. This approach generally re- 
stricts itself to a station analysis of the fuselage. 

The vehicle synthesis and optimization procedures used in 
the ACSYNT program have evolved using empirically based equations. 
Because the user typically does not have the expertise to gener- 
ate all inputs alone it has been decided that the more detailed 
structural design programs would not be executed from within the 
synthesis loop. The more detailed structural designs would be 
made outside of the vehicle synthesis loop either to calibrate a 
point design or to generate analytical data which may be used to 
modify the statistically derived equations. 
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During this contract period, a program for aircraft struc- 
tural synthesis (PASS) was developed by Erwin Johnson under the 
sponsorship of the National Research Council. This beam finite- 
element program has initially been written to analyze the aero- 
elastic effects of wing flutter and redesign. It is being devel- 
oped with the full vehicle analysis capability as a goal. 

The approach taken in developing this fuselage design code 
incorporated each of these items. Rather than restrict the ana- 
lytical procedures to simplistic approaches, the more detailed 
structural analysis techniques (Level II) employed in the SWEEP 
program were adapted for the estimation of basic fuselage shell 
weights. The program structure was then written to accept the 
shear, bending moments and external loads from the PASS program, 
or as direct input. It is anticipated that this approach will 
provide more detail and potential in which to assess advanced 
technology. This should satisfy the third function of providing 
analytically derived data which may then be used to check point 
designs or modify existing statistical equations. 

In developing the computer code, certain portions of the 
fuselage design code were extracted directly from the SWEEP pro- 
gram. The modif ications made to the original code were to improve 
the readability of the code and to incorporate the geometric des- 
criptions used by the ACSYNT program. In certain instances 
failure criteria, such as sizing the fuselage for acoustic fa- 
tigue, have been deleted because of lack of input from the ACSYNT 
program. Their deletion from the code will be marked for future 
reference. 

The details of theoretical background behind the develop- 
ment of the computer are described here. The basic flow of the 
program and the equations used to size the primary structure are 
also outlined. An example case for the C-141 follows. 
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Approach to Weight Estimation 

The basic approach taken in computing fuselage weights will 
be to compute the items listed in figure 14. The weights pre- 
dicted here are categorized as basic or secondary structure ac- 
cording to the definitions of MIL-STD-1374 and are listed by 
line items on form AN-9102-D. This is the basic military stan- 
dard form of weights classification. The civil transport aircraft 
weight reporting format essentially follows this form with minor 
modifications for special items. 

All items in this weight statement are computed either by 
direct analysis or by empirical equations. Most analytical 
equations used to predict component weights with the exception 
of the sizing of the major frames have been excerpted from ref- 
erence 2. The major frames computations have been left out until 
such time as integrated approach is implemented which incorpor- 
ates the interaction between fuselage, wing, tail, and landing 
gear loads. They will be accounted for by the use of empirical 
equations . 

Of the items listed under BASIC STRUCTURE the weights of 
bulkheads, minor frames, covering, stiffeners, and longerons are 
computed analytically by estimating the required material distri- 
bution at up to twenty synthesis cuts. The remaining items and 
all SECONDARY STRUCTURE items are computed individually from 
empirical equations. Three methods are available for computing 
most secondary structural weights. They may be input, estimated 
from rule-of-thumb weights, or calculated from geometric infor- 
mation for fighters, transports, and bombers. Most calculated 
weights are estimated based on a unit weight of the items. 

Weight of the primary structure is estimated by computing 
the required material gages at each synthesis cut. An outline 
of procedure used to generate these estimates follows: 


1. Input geometry and loading conditions. 

For each synthesis cut: 

2. Determine maximum loading condition and material prop- 

ties. 

3. For a given stringer and minor frame spacing compute 
unit weight of minor frames, cover panels, and longerons or 
stringers . 

4. Increment the stringer or minor frame spacing if desired 
to search for minimum weight spacing. Return to (3) until mini- 
mum unit weight is obtained. 

5- Compute minimum weight pressure bulkheads if required. 

At completion of synthesis cuts: 

6. Summarize primary structural weights and apply non- 
optimum weight factors. 

7. Compute secondary structural weights. 

8. Print summary. 

FUSELAGE GEOMETRY DESRIPTION 

All geometric dimensions used in the analytical equations 
are obtained from the description of the external geometry. The 
basic technique of inputing the geometric shape employs the use 
of a restricted set of automatic geometry descriptors (AGD's) to 
describe the basic shell size and shape. The descriptors derived 
here are restricted to circular and double-lobed circular shapes. 
Reference 2, which may also be implemented, employs a basic 
rounded- rectangle as the primary AGD. 

The external geometry is defined according to the control 
option, IGM. For IGM = 1, the fuselage is defined according to 
the ACSYNT program vehicle synthesis AGD to have a circular cross 
section and a Sears-Haack area distribution. For IGM > 1, the 
particular AGD description of the local cross section is input 
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at up to ten body stations. Of these geometry cuts, the first 
and last define the nose and tail. Intermediate cuts may cur- 
rently be defined at up to eight stations. Sharp geometric 
changes, such as occur forward and aft of duct inlets, are des- 
cribed by double cuts immediately forward and aft of the shape 
transition. The standard body axis system with the x-axis posi- 
tive out the tail and the z-axis positive up is used to define 
the vehicle coordinate system. 

Structural sizing is performed at up to 19 synthesis cuts 
(XO) . These cuts are located at NC stations along the longitu- 
dinal axis of the vehicle between the nose and the tail. These 
locations do not have to correspond to the locations defining 
the external geometry. However, they must be used to define the 
boundaries of such structural features as cutouts, bulkheads, 
and locations of particular interest to the user. Shell dimen- 
sions at these synthesis cuts are obtained by linear interpola- 
tion between AGD descriptions of the external geometry. A des- 
cription of the parameters that make up the circular and double- 
lobe circular AGD's follows. For a description of the rounded- 
rectangle see reference 2. 

In the discussions that follows, the term "cut" refers to 
a synthesis cut at which the various structural details such as 
cross-sectional geometry, loads, and material gages are evalua- 
ted. Similarly, the term "segment" refers to vehicle properties 
computed between two synthesis cuts. The following section des- 
cribes the geometric properties and equations derived for the 
Sears-Haack body input (IGM = 1) and the circular and double- 
lobe circular shapes (IGM = 2 ) . 

The Sears-Haack body is assumed to be circular in shape and 
is partitioned into three sections: a nose, a constant section 

body, and a tail. The input parameters (IGM = 1) defining the 
body radius and length at any cut are: the body length (BODL) , 


the body diameter (BDMAX) and the nose (FRN) and tail (FRAB) 
fineness ratios. Figure 15 is a picture of the Sears-Haack body 
and circular cross section. The radius of the cross section at 
any station is defined by: 


r = 


BDMAX 


3/1 


[iliy-SL] , for 0< 


- BDMAX - 


< FRN 


(40) 


where: 


5 “ 2 FRN BDMAX 


BDMAX 

2 


for 


FRN < 


— - — < 1 
BDMAX x 


FRAB 


(41) 


r = 


BDMAX f £(1 ~ £) 


3/4 


for 1 - FRAB < 


- BDMAX - 


(42) 


where : 


BODL - x 
2 FRAB BDMAX 


The floor location in this description has been specified 
at a uniform height above the bottom of the cross section of 
maximum diameter. It is input as the ratio (ZFLRU) of the floor 
height to maximum body (BDMAX) diameter. The locations of syn- 
thesis cuts (XO) are input as the distance from the tip of the 
nose. All remaining section properties are either input for the 
synthesis cut location or are interpolated using the description 
of the body radius. 

The input parameters that define the double-lobe circular 
cross section in figure 16 ( IGM = 2) are the scaling factors 
(RFLB and SCALEB) , the perimeter correction factor (PERI) , the 
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nondimensional values of the upper lobe radius (Al) , the lower 
lobe radius (A2) , the vertical offset of the centers of the upper 
radius (El) , and the lower radius (E2) . These values are non- 
dimen sionali zed according to the scaling factors 


Al 


Rl 


SC ALE B] 
RFLB 


(43) 


The input cross sections are defined at up to NGM (less than or 
equal to ten) stations (XI) . The location of the floor (ZFLR) 
is either input as the ratio of floor height above the bottom to 
the local body diameter for circular sections or computed as the 
height at the intersection between the upper and lower lobes for 
noncircular sections. The perimeter may be adjusted by multi- 
plying by the correction factor, PERI. The locations of synthe- 
sis cuts (XO) are input with respect to the user defined vehicle 
reference axes. The double-lobe cross sectional description re- 
duces identically to the circular section for either El s E2 
or for A2 = 0. In the latter case, A2 is set equal to Al 
and E2 is equal to El. 


SECTION GEOMETRY 

The structural sizing techniques in this program follow the 
methods employed in reference 3. For structural sizing the cross 
section is partitioned into four shell sectors representing the 
upper, lower, and two sides. For the purpose of simplicity, the 
shell sectors are assumed to be symmetric about the centerline 
of the aircraft. The section quantities estimated for the struc- 
tural geometry in routine FSECTN are: the radii of the sector 

(R cu , R cl * and R cs ) and the sector arc lengths BU, BL, and BS) 
for the upper, lower and two sides, the total perimeters (PER), 
the area of the cross section (ARCS) , the maximum section depth 
(DF) , the maximum section width (WF) , the ratio of the floor 
height to section depth (DKHT) , and the width of the floor at 
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the section (W pLR ) . Figure 17 shows the relationships among 
these various geometric quantities for the double-lobed section. 

For program efficiency the upper and lower lobes are des- 
cribed internally by the angles from the vertical to the deck 
intersection. The upper lobe is defined by the arc swept out 
by ALPU measured clockwise from the positive centerline. The 
lower lobe is defined by the angle ALPL measured counterclock- 
wise from the negative centerline. The structural sectors are 
similarly defined by the angles 0 and 6 0 . They are limited 
to 45° of the upper and lower lobe angles. If the section is 
circular the sectors are still defined by the location of the 
floor. If no floor exists, the upper lobe is defined by 
ALPU = 3tt/4 and the lower lobe is set to ALPL = tt/4. The equa- 
tions used to generate these parameters follow. 

The following quantities are defined for the double-lobed 
cross section. The maximum section depth is: 

OF = R cu + R cl + EU - EL (44) 

The deck height to section depth ratio is: 

R - R cos 0 + EU - EL 
DKHT = — == (45) 


where : 


cos 


[(EU - EL) 2 + R cu 2 - R cl 2 ] 
2 R CU (EU - EL) 


(46) 


The floor width is: 

WFLR = 2 V r cL 2 " ( DKHT#DF “ R cl )Z (47) 

The angle of the upper lobe arc measured from the positive 
centerline is: 
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ALPU 



sin 



- R cl + DKHT • DF 
R CU 


(48) 


The angle of the lower lobe arc measured from the negative 
centerline is: 


ALPL - -• 


- 1 

sin 



- DKHT»DFl 


CL 


(49) 


The maximum section width is: 

f R cu »sin £ min (ALPU,j)j 

WF = max l 1 (50) 

I R CL »sin j min (ALPL,j)J 

If the section is circular and no floor exists the geometric 
properties may then be defined as: 


DF - 2 • R 


CU 


DKHT = 0. 


WFLR 

ALPU 

ALPL 

WF 


= 0 . 

3tt 

4 

IT 

“ 4 
= 2 • R 


CU 


(51) 

(52) 

(53) 

(54) 

(55) 

(56) 


The section sectors used in syntehsis are similarly defined 
by the angles 0 and 6 C measured from the centerline to the 
boundary. These sectors partition the perimeter into three 
regions according to the types of loading they will be carrying. 
The upper and lower sectors are designed for fuselage bending 
loads and stiffness while the two side sectors provide shear 
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is measured from 


strength and lateral stiffness. The angle 6 u 

the positive centerline to the sector boundary. The angle 9 £ 

is measured from the negative centerline to the sector boundary. 

The angles 6 and 6. are defined as: 
u t 

9 U * min (ALPU, £) , and (57) 


0 £ - min (ALPL , J) (58) 

The choice of upper bound of n/4 is arbitrary. For 
circular sections this partitions the perimeter into four equal 
sectors . 

The fuselage perimeter (PER) and sector arc lengths (BU, 

BL , BS) are then defined as: 


PER = 

2(R CU -ALPU 

+ R cl *ALPL)*KC 

(59) 

BU = 

2 "W 9 u 


(60) 

BL = 

2 - r cl \ 


(61) 

BS = 

i(PER - BU 

- BL) 

(62) 


The radius of curvature of the side is approximated from 
the weighted average of radii of upper and lower sector multi- 
plied by their fraction of the side sector perimeter. This is 
used to provide a measure of effective radius when estimating 
the curvature corrections in various sizing criteria. The equa- 
tion for the effective side radius is: 


R 


CS 


-R CU 2 (ALPU - 0 u ) + R cl 2 (ALPL - 0 £ ) " 

— 


(63) 


The centroid of the cross section is approximated as the 
centroid of a shell perimeter of unit thickness. The lateral 
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centroid is assumed to lie on the centerline. The vertical cen- 
troid which defines the neutral surface is: 


ZO - 


2 R CU (EU«ALPU + R cu *Sin ALPU) + 2R £L (EL* ALPL - R CL *sin ALPL) 

PER 


(64) 


The location of the center of curvature for the upper and lower 
lobes relative to the vertical centroid is: 

Zy - EU - ZO (65) 

Z. » EL - ZO (66) 


The cross sectional area is: 

ARCS = R cu 2 *ALPU + R cl ? ‘ALPL + j W FLR (EU - EL) (67) 


Segment geometry is defined for four types of segments: 
the nose, the tail, the normal segment, and the sharp transition- 
al segment. The four segment properties defined in routine FSECTN 
are the length (DELX) , the surface area (SF) , the internal volume 
(VOL) , and the center of mass of the volume (XBAR) . The segment 
properties of the nose are: 

DELX = XO, - XI j (68) 


(2r + r ) 

XBAR = XOj - DELXj V y y (69) 

SF = ttUj + r 2 ) *V DELX 2 + (r t - r 2 ) 2 

(70) 

VOL * J DELX ( r 2 + r 2 rj + r 2 ) (71) 
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where 


r t = radius at nose 

PERj 

r = - - . 

2 2" 

The segment length at the tail is: 

DELX = XI NGM - X0 NC 

The equations for XBA SF, and VOL at the tail are the same 
for the nose with the exceptions that: 


r = 


per nc 


i 2ir 


and 


r„ = radius at tail 


The properties for intermdeiate segments are: 


DELX. = XO. - XO. , 

3 3 3 ~ 1 

XBAR . = i (XO. + XO. ) 

3 2 3 3-1 

SF. = tt ( r + rjVoELX 2 + (r, - r ,) 2 

I 2 w 2 1 

DELX . , 

VOL . = ( ARCS . + ARCS . , + V ARCS . ARCS . . 

3 3 \ 3 D “ 1 v 3 D" 1 


where 


r i = 


PER. 
IZl 

2 TT 






(72) 

(73) 

(74) 
as 

(75) 

(76) 

(77) 

(78) 

(79) 

(80) 

(81) 
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and 


(82) 


r 


2 



If the segment is a sharp transition (DELX 2 inches) the sur- 
face area and volume are approximated by: 

SF . = PER. DELX . (83) 

3 3 3 

VOL. = DELX. ARCS. (84) 

3 3 3 

The unit inertias for the fuselage and its contents are 
generated in routine INERT. The unit inertias of a segment are 
the inertia per unit weight of the fuselage shell and its con- 
tents about its centroid. The inertia is computed assuming the 
segment is a solid cylinder with the average properties of its 
ends. The unit inertia is estimated by dividing the local volume 
inertia for a solid cylinder of average dimensions by its cross 
sectional area. The pitch (UIY) , roll (UIX) and yaw (UIZ) unit 
inertias are calculated by assuming: 

1. The weight of the fuselage and its contents are uniform- 
ly distributed within the volume of the segment. 

2. The center of mass is at the centroid of the segment. 

The inertias per pound of weight for the double-lobed cir- 
cular shape are defined for the normal geometry transition from 
the following: 



( ALPU . + ALPU . , 

3 3 " 1 


(ALPL . + ALPL . ) 
3 3 


R CIL + R CUj_j 


(85) 

( 86 ) 
(87) 
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UIY . 
3 


= ( 2r' ± Z, + i r' Ma + ^ r’z'-sin a 

1 u M2 u 8uju 3uu u 


+ r ' sin 2a 
16 a u 


+ 2r' n 


1 7 ' 2 + I r ' 2 

2 Z Z 8 l 


a , 


+ 3 r £ Z £ r,in a £ + 16 r £ sin 2 


a *lM 


DELX‘ 

12 


UIZ . = 
3 


1 . 4 

4 r u 


a u - j sin 2a u 


. 1 > 4 ( 

+ ? M 


“t - 2 sin 2a l 


•fl , DELX. 

A" + 12 " 

3 


DELX 2 

UIX. = UIY. + UIZ. - — 7— L 
3 3 3 6 


For sharp transition segments the unit inertias are based 
on the properties at the aft end of the segment. 

For the nose and tail segments, the equivalent section 
radius is used to calculate the inertia. The nose segment in- 
ertia is defined as follows: 




9 


and 


(96) 


UIX, 


3(r? 

10 (r ’ 



- rj) 


UIX » 27 

UIY i = UIZj = — (XO, - XBARj) 


(97) 


The inertia of the fuselage is computed by summing the 
segment inertias about the aircraft center of gravity. 

TORSIONAL GEOMETRY 

The internal geometric arrangement is required in order to 
determine the thickness to satisfy any torsional rigidity (GJ) 
requirements at a given synthesis cut. Evaluation of the tor- 
sional capability of the shell is based on the presence of a 
closed torque cell. The presence of an open torque cell is 
considered structurally inefficient, and will not be considered 
in these calculations. The values calculated here are also used 
in the evaluation of pressure bulkheads. 

The thickness required to satisfy a given torsional rigidity 
is derived from the torsional constant for thin-walled closed 
sections of general shape. The torsional constant, J, is defined 
as follows: 



ds 

t 


The thickness required is estimated for a uniform thickness (TGJ ) 
around the perimeter of the torque cell as: 


TGJ 


(GJ) 


reqd 

4A 2 G 


PER 


(99) 
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where: 


A = enclosed cross-section area 
D = panel shear modulus 

PER = peripheral length of enclosed torque cell 

The only geometric parameters required in the evaluation of 
the required thickness are the enclosed torque cell peripheral 
length (PER) and cross sectional area (A) . Since the presence 
of cutouts is assessed, the internal geometry must be capable of 
establishing a consistent torque cell. The primary internal 
partition that may be used to define the cell configuration is 
the horizontal deck. The presence of a cutout defines the ori- 
entation of the torque cell, upper or lower. The horizontal 
deck, then, is the structural member that establishes the closed 
torque cell shape. All decks and upper and lower cutouts are 
bounded by synthesis cuts. Therefore, at any synthesis cut, 
there may be two different conditions defining the structure 
forward and aft of the cut and two corresponding required thick- 
ness values. These differences are evaluated by calculating 
torque cell data on both sides of each cut. 

Certain geometric combinations are not compatible with the 
assumption of a single closed section. The program approach for 
these arrangements is as follows: 

1. Should decks exist without any cutouts, the external 
section geometry is used to define the torque cell. The influ- 
ence of the deck upon torsional stiffness will be ignored. 

2. Should cutouts exist without any decks, the external 
section geometry is used to define the torque cell, and the loss 
of torsional stiffness due to cutouts will be ignored. 

3. Should both upper and lower panel cutouts exist in the 
presence of a deck, the section above the deck is used to de- 
fine the torque cell, and the loss of torsional stiffness due to 
the upper panel cutout will be ignored. 


67 


Figure 18 depicts the different internal arrangements and 
the corresponding geometric variables immediately forward of a 
cut. Variables in the calculations are: 


ACRS 

PER 

ANTF, 

ANTA 

Total shell cross-section area at cut 

Shell external perimeter at cut 

Torque cell cross-section area immediately 

PERF, 

PERA 

forward and aft of the cut, respectively 
Torque cell peripheral length immediately 

PRDF, 

PRDA 

forward and aft of the cut 
Deck peripheral length immediately forward 

DEPF, 

DEPA 

and aft of the cut 

Depth of the torque cell immediately forward 

WIDF , 

WIDA 

and aft of the cut 

Width of cell between two walls (beaming 


distance) immediately forward and aft of the 
cut 


For the presence of a pressure bulkhead at a synthesis cut, 
the above variables define the geometry of the pressurized 
compartment. In these instances the width, depth, perimeter, 
and area describe the dimensions of the bulkhead. 

In addition to tne basic closed cross section, two arrange- 
ments that incorporate cutouts in the upper and lower covers are 
defined. The equations defining the structural geometry of the 
torque box for a cutout in th i lower sector shown in figure 19 
are : 

PERF = 2R CU *ALPU + W FLR (100) 

DEPF = DF ( 1 - DKHT) (101) 


f ^ * R cu 

WIDF = l 

, if 

ALPU > 

90 

l 2»R_ *sin ALPU 

, if 

ALPU < 

90 

v CU 





( 102 ) 

ANTF = ALPU*R CU 2 - | W FLR < R CU " DEPF) (103) 
PRDF = W FLR (104) 
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Similarly, the equations defining the structural geometry of the 
torque box for a cutout in the upper sector shown in figure 20 
are : 


PERF 

DEPF 

WIDF 

ANTF 

PROF 


= 2R CU *ALPL + W flr 




(105) 

= DF-DKHT 





(106) 

r 2R c L -Sin 

ALPL , 

if 

ALPL < 

VO 

O 

O 


"'- 2r cl 

9 

if 

ALPL > 

90° 

(107) 

- ALR 1- R CL 

I W FLR 

(r cl 

- DEPF) 


(108) 


‘ W FLR (109) 


The minimum material thickness required to satisfy the 
torsional stiffness requirements (GJRD) may then be written for 
the forward segment as: 


GJRD » PERF 
4ANTF 2 *G 


( 110 ) 


The properties for the aft segment are defined in the same 
manner . 


STRUCTURAL GEOMETRY AND INERTIAS 

Determining the precise internal loads distribution is not 
within the scope of this program. The maximum bending and shear 
stresses in the cross section are estimated from an equivalent 
beam analogy. If certain assumptions are made, the internal 
load distribution can be approximated solely on the basis of 
structural geometry. The primary assumptions are first that the 
centroid of the material area moment lies at the centroid of the 
unit cover; second, that all sizing elements within a sector are 
the same; and third, that the bending contribution of the skin 
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in a postbuckled state is neglible compared to the contributions 
from the longerons or stringers. 

From beam theory, the maximum shear flow occurs along the 
axis of the centroid and can be expressed as: 


q 



( 111 ) 


where : 


V = vertical shear at the synthesis cut 
q = shear flow in pounds per inch 

Q = area moment of the axial elements in one quadrant 
I = total area moment of inertia of all the bending elements 

Similarly, the shell bending stress at any vertical coor- 
dinate relative to the centroidal axis is: 


a 


M (Z - ZO) 
I 


( 112 ) 


where : 

M = bending moment at the synthesis cut 
Z = vertical coordinate of the member 
I = total area moment of inertia 
ZO = vertical location of the area centroid 


The maximum positive and negative stresses occur for the corres- 
ponding external distances from the centroid. 

The area moment, Q, is defined as the summation of the 
axial elements in the upper quadrant: 


Q = 


Zda = A LfJ lz + A Tfi £Z + A tt lZ 


LS 


IT 


+ T cu ! zda + T cs 


Zda 


( 113 ) 


where: 
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A 

A 

A 

T 

T 


LU 

LS 

IT 

CU 

CS 

ds 


Iz 

The 


I 


■ area of each upper longeron or stringer 

• area of each side stringer 

* area of each secondary longitudinal member 
= thickness of upper sector cover 

= thickness of side sector cover 
= incremental panel length 

= summation of vertical coordinates of elements 
area moment of inertia is: 


Z 2 dA = A lu IZ* + A ts [Z 2 + A tt [Z 2 + A tt lZ 


LS 


IT 


LL 


’cu z ' ds + T CS z2ds + t cl z2<3s 


(114) 


where : 


A ll = area of each lower longeron or stringer 
T^ t = thickness of lower sector cover 


The approach taken in this program is to compute the area 
moments and inertias for a unit thickness or area of the con- 
tributing elements in each sector. This contribution is a 
function only of the arrangement of the bending elements: covers 
longerons, and stringers. The initial approximation to the 
maximum shear is based on the geometry. 

The arrangement of structural members for the double-lobed 
circular sector is shown in figure 21. Two optional structural 
arrangements are depicted for longeron and stringer construction 
Longeron construction is characterized by four primary longitu- 
dinal members at the boundaries of the sectors which sustain the 
bending load. Stringer construction, on the other hand, is 
characterized by longitudinal stiffeners at an even spacing, 
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BSTR , around the perimeter of the shell. The unit inertias for 
each sizing member in each sector are calculated in subroutine 
IL0NG1. For the purposes of consistent notation the symbol I y 
will refer to vertical inertia about an axis parallel to the 
y-axis. Likewise, I z will refer to the lateral inertia about 
the z-axis. 


The equations used to calculate the ur.it vertical bending 
inertias about ZO for the upper, lower, and side sectors are: 


fix.1 

= 2R_.. 

* / 

t T cJ 

CU 

l 


7 2 . 1_ p 2 

Z U + 2 R CU 


'CL 


= 2R 


CL 


[(v * i 


R 


CL 


6 U + 2 V R cu sin °U * 7 R cu !sin 2 e i 
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e i * 2 V r cl sin e i + 7 R ci> in 2 e t 


(116) 
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+ 2r cl | 
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^ T ci> 
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^ t clJ 


ALPU+2Z U *R CU sin ALPU + i R cu 2 sin 2ALPU 


ALPL + 2Z t *R_. sin ALPL + t- R 2 sin 2ALPL 
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The equations used to calculate the unit lateral bending 
inertias about the centerline for the upper, lower, and side 
sectors are: 


1 z 1 
• T CIP 

= R cu 3 

CD 

c 

1 

2 
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2 0 

U J 
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( 119 ) 
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R 3 

R cu 


ALPU 


1 

1 


sin 2 • ALPU 



r cl 3 


ALPL - * Sin 2 • ALPL 


) - (t~) 

> ^CL' 


( 120 ) 


The distance of the extreme fibers of the upper and lower 
covers from the neutral axis are defined for vertical bending 
as follows: 


Z = Z„ + R_ t , 
max U CU 

Z . = Z T - R„ t 

min L CL 


and 


In the presence of cutouts, a cutout longeron is positioned 
along each side of the hole to carry the bending load around the 
removed structure. The cutout longerons are assumed symmetric 
about the vertical centerline and computed wherever the effective 
cutout widths RTU and RTL for the upper and lower sectors are 
nonzero. The unit inertias of these longitudinal members are 
proportional to the square of their distances from the neutral 
axes. The horizontal and vertical distances of the upper cut- 
out longeron ( Y cu » z cu^ are as: 


Y C u = R cu‘ sin 6 ' and 


( 121 ) 


z cu = z u + R cu cos 0 


( 122 ) 


where 6 = 1/2 (^u^CL^ an 9 u l ar location of the longe- 

rons measured from the vertical centerline. The horizontal 
and vertical distances from the neutral axes for the lower 
cutout longeron (Y_ T , Z_ T ) are defined as: 

Li.) Lij 


73 


Y CL ■ R CL ,Sin 9 
Z CL ■ Z L * R CL COS 9 


( 123 ) 

( 124 ) 


where 6 * 1/2 (R^/R ) measurec * ^ rom the negative centerline, 

The unit inertias for longitudinal stiffening elements are 
computed for two structural arrangements: stringers and longe- 

ron construction. The stringer unit inertias are estimated 
assuming the stringer spacing, Bg TR , is sufficiently small in 
comparison to the perimeter that the effective unit inertias can 
be distributed as an equivalent thickness. The unit inertias of 
the stringers may then be written in terms of the unit inertias 
of the cover as : 

Upper, Lower and Side Vertical unit inertias - 



( 125 ) 


( 126 ) 


( 127 ) 


( 128 ) 
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The maximum shear stress due to bending is computed from 
the maximum shear flow at the neutral axis. The area moment , Q, 
used to calculate the shear flow may be estimated from the geome- 
try of only one quadrant of the shell. The equations for the 
vertical and lateral unit area moments for the upper lobe string- 
ers in figure 21 are: 


* R 


(2^*0 + R cu sin 0) 


LU ; 


V 


cu 


B 


STR 


( 131 ) 


LU 


= P. 


2(1 - COS 0 ) 


CU 


B 


STR 


( 132 ) 


where 9 is the angle from the centerline to the intersection 
of the neutral axis, ZO, with the upper lobe. If 0 is greater 
than ALPU, the lower lobe geometry is used. 

The unit inertias for longeron construction are approximated 


for primary and secondary longerons as: 




where 


Z LNG Z U + R CU COS 0 LNG 


2 (Z L “ R CU* COS 6 LNG ) 


( 135 ) 


SEC 2 


(136) 


where is the angular location of the primary longerons 

measured from the sector centerline. The vertical location may 
optionally be located by the ratio of the vertical height to 
total fuselage depth. N LNG anc * N SEC are num ^ ers ° f primary 

and (if requested) secondary longerons. Secondary longerons are 
located at half the height of primary longerons. Similarly the 
unit lateral inertias ares 


r*- ‘T 


y2 \ (R CU Sin °LNG ) 


^ (r cl sin 9 lng ) 2 


i[ (R C u Sin e SEC )2 ] 


( 137 ) 
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Only the primary longerons are assumed to contribute to the 
shear flow. The equations for the vertical and lateral unit area 
moments for the upper lobe longerons are defined as follows: 


(JL| = 

a t „ 

V T.Il/ 


LU' V 


'LNG 


( 138 ) 


'LU J 


= R cu sin 6 lNG 


J 


The properties of the internal and structural geometry for 
the rounded rectangular shaped cross section are given in refer- 
ence 2 . 


Shell Cover Design Criteria 

During the sizing of a given shell segment a number of 
assumptions are made to simplify the analysis. For the synthesis 
of shell structures in bending the circumference of the segment 
is assumed to be partitioned into four quadrants: two sides, 

and upper and lower sectors. All bending loads are assumed to 
be carried by the upper and lower sectors, while the sides are 
assumed to carry all shear loads. The contribution of the long- 
itudinal members within all sectors is included when determining 
the net bending moment carried by the section, and when deter- 
mining the maximum shear stress along the sides. Frames are 
sized to provide shell stability, and stringers are sized to 
carry both bending load and provide buckling stability. The 
basic covers are also checked for pressure integrity. The sizing 
criteria examined are broken into four types: shear criteria, 

bending criteria, stability criteria, and pressure design cri- 
teria. The following constraints will be used to size the fuse- 
lage cover elements: 

Cover Element Sizing Criteria - 

(1) Minimum gaqe - T_ > TC . (139) 

l — mm 
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(2) Ultimate shear stress - T_ >_ p — TF ' 

* su k 

(3) Shear buckling - f < F 

5 “ S C i 

(4) Postbuckling - f < F_ 

s “ s allow 

= F + sin a cos a(C_ F. 
scr R tu 

a = 45° 


(140) 

(141) 


F 

scr 


(142) 


(5) Cabin pressure 


(M.S. 


(6) Diaphragm pressure 


PR 

T c > _= (M.S.) 

2.0 manned 
1 . 5 unmanned 


1.3769 b 
> 


2.484 

P 

4.467 

o 


(143) 

1.984 

(M.S. ) 
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1.646 


b 


0 897 0.394 

P E C 

1.288 

o 


(m.s . ; 


(144) 


(7) Local panel flutter (M>1) 



T F = 


B 


F (M) • E 


1/3 


'BNDRY 


(145) 


;p n = 0.5551841-0.1686944 


f l! 
w 


+0.02169992 


+0.00096394 = 


L' 1 


wj 


( 146 ) 


where 


f - sector shear stress = q/T 

S V" 

T l - panel thickness at landings 

q - maximum shear flow at neutral axis; q = V^Q/1^ 

F - ultimate shear stress 

su 

F - critical buckling stress 

scr 

F. - ultimate tensile stress 

tu 

R c - radius of curvature of shell sector 

o - material allowable stress 

p - shell limit pressure, psi 

q^ - cover flutter critical Mach number | 

F(M) - Mach number correction obtained from reference 2, 

Vol VII, figure 14. j 

C BNDRY ~ boundary layer correction j 

Only the side sector cover panels are sized for the maximum ] 

shear stress. Therefore, only the side panel thicknesses are 
checked against the ultimate shear stress, shear buckling, and 
postbuckling requirements. 

MINOR FRAME SIZING CRITERIA 

During the synthesis of a given shell cut, miner frame or 
ring design is carried out at fixed frame spacings (SFRM) and j 

the total shell weight determined. A search is then made for 
the spacing which predicts the least weight. Figure 22 is a 
drawing of the assumed geometric model used in the analysis of 
the minor frames. Several simplifying assumptions have been 
made in the geometry to minimize the number of design variables. 

The frame depth (FD) and cap width (BFCM) are input, by the user 
or from default values. Only the ring thickness (t r ) is sized. 


| 

j 

i 

I 

*3 

i 

j 
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The geometric properties of the ring may then be written 
for the area, area moment of inertia, and radius of gyration as: 


Ring area: 


♦ rp] 


A = t 4 BFCM + 

r r [ 2 J 


(147) 


Ring moment of inertia: 

I = t [bFCM FD 2 + \ BFCM 3 - BFCM 2 FD + ?? 
xx r 3 24 


(148) 


Radius of gyration: 


(149) 


The equations used to size the ring thickness are explained 
in reference 2. A summary of the equations and constraints used 
in subroutine MINFR is given here. The constraints used to 
design the ring thickness are: 

1. Minimum gage requirement (TFCM) : 


t > TFCM 
r - 


(150 


Shanley's general shell stability requirement 


. max 1 tt | CF 

r - (I 

E SFRM -~ 

r l t r J 
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= magnitude of maximum bending moment 


effective body diameter at the synthesis cut 

Shanley's coefficient to prevent general shell 
instability (= 1/16000) 


E f = ring modulus of elasticity 

SPRM = minor frame spacing 

J xx 

— = moment of inertia per unit thickness 


3. Forced crippling requirement for shear panels: This l 

constraint is checked when the cover shear panel is considered 

I 

to be postbuckled (f F ) . The forced crippling requirement j 

is obtained iteratively between satisfaction of equilibrium of 
the net loads carried by stringers, rings, and skins and the 
compatibility of stresses and strains. This is contingent upon j 

estimating the diagonal tension angle, a, at which the load is j 

effectively carried through the skin. j 

The diagonal tension factor for curved plates is: ] 


where 


K 


tanh 


0.5 + 300 


TC • RLD 


CS 


log 


10 


scr 


R CS 

RLD 


cover thickness 

shell radius of shear panel 

SFRM BP AN ' 

maX [BPAN * SFRM' 


(152) 
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BP AN 


stringer 

longeron spacing 


Where stringer and longeron equation differences occur, the two 
will be presented side by side. The initial estimate for diag- 
onal tension angle is: 


Stringer 


Longeron 



A = 0 (153) 


Since the stringer/longeron area is not yet known, an initial 
estimate of the required area to satisfy minimum area, ASTRM, 
and/or maximum bending load is made. 

Stringer Longeron 

A > ASTRM A > ASTRM (154) 

s — s — 


M ^ 
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l. 

cy 


(155) 


Then 


a PDT 


5 - + 0.1443 A 

4 

1 + 0.175622 A + 0.013411 A 2 


(156) 


And 


a 


a PDT K 


2 5 


(157) 


The load in the stringers is distributed evenly between stringers. 
The load in the longerons i s carried only between the upper and 
lower longerons. The maximum average stress may then be written 
for each as. 
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The maximum induced stress in the ring versus the maximum 
average stress is defined in terms of the panel aspect ratio as 


RG MAX _ 


1.0 


if §FRM > 1 

i ± nnim x ^ 


RG 


1.0 + 0.78 ( 1-K ) -0 . 65 ||^ (1-K), if 


BPAN 

SFRM 


BPAN 


BPAN - 1,2 
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The allowable ring stress is 


RG 


= A K 


V' 


allow 


TC 


1/3 


(160 


where 


r 


A = < 


(0.18695 + 0.00075238 R cs ) 


.E . 
v c ; 


1/9 


, 0 R cs < 151.586 


10. 30100 


rE ) 1/9 

r 


R cs > 151.586 


(161) 


and 


G = 




5.88 


+ 0.002 

E r 


(162) 


The solution for a and t is iterative. The following 


calculations for a are carried out three times to obtain an 
approximate value. 
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Strain in covers 


= =~ — ■ % — + sin 2 a (1 - 


'c E sin 2 a 
c 


K) (1 + y) 


( 163 ) 


Strain in ring: 


e = 

r 


RG 

E 


( 164 ) 


Strain in stiffener: 


r _ f ST 
E S E 

s 


( 165 ) 


the value of a is computed from: 
Stringer 


Longeron 


tan a 
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tan a = 


e - e + 
c r 


e - e 

c s 


1 

Cbpan' 

2 

e . e + 1 

sfrm' 

24 

i* R CS ' 

1 

c e r 8 

■ R CS > 


tan 


(166) 


and 


a = tan 1 a 

The iteration on the calculation of a, either returns to equa- 
tion 163 or terminates after three cycles. 

The ring thickness required is computed from the following 
quartic equation: 

Wb + X c> 3 - X a = 0 <167) 

Newton's method is used to solve for the ring thickness. If 
t is greater than the value estimated for stability and minimum 
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gage, the procedure starting with the a iteration is repeated 
a second time to update a. The coefficients X^, and X c 

are: 


X 


a 


ivf tan ai 

(-SH 


RG MAX 

T 


RG 


)] 


3 


K 



FD 

4 BFCM + 


1 + 

’ FD ' 

i 2 v 

1 

J 

SFRM T c 


X c = 0.5(1 - K)(j£] 


(168) 


(169) 


(170) 


where 



D 


stringer/longeron modulus of elasticit' 
cover modulus of elasticity 
ring compressive yield strength 


moment of inertia contributions per area of 
upper and lower longeron or stringers 

total body depth at cut 


SHELL BENDING CRITERIA 

The lv. gitudinal material required to resist bending is 
determined from the sum of the components that may carry bending 
loads. The covers are first sized to satisfy pressure require- 
ments. The bending contributions of the cover material are then 
determined for tensile and compressive sides. On the tensile 
side, the entire cover is considered effective. On the compres- 
sive side, reductions are made for the effective width of the 
skin carrying moment. The difference between the total moment 
required and the moment carried by the cover is used to size the 
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longerons or stringers. The contribution of the side covers to 
resist bending loads .is neglected. The presence of cutouts is 
predicted by elimination of cross sectional material where indi- 
cated, and by material addition in other areas. Only the maxi- 
mum up-bending and down-bending loading conditions are used to 
size the material. In ac'dition, bending stiffness (EI^ anc * EI Z ) 
requirements are checked. 

Figure 23 is a sketch of the geometric configuration of 
the assumed model of the longitudinal member used in evaluation 
of the forced crippling strength of the shell. Several simpli- 
fying assumptions have been made in the geometric proportions 
to reduce the number of variables. Stiffener (longeron) geome- 
try is set up with flange width equal to web height, HSTR. The 
stringer flange to height ratio, RSTRII, and web height have 
programmed default values which may be overridden by user input. 
The area, first moment, and the moment of inertia of the stif- 
fener are: 

Are-i for forced crippling: 


AFC = tgHSTR ( 3 + RSTRH) 


( 171 : 


First moment and centroid: 


AFCY •= t HSTR (0.5 + RSTRH) 
s 


e = HSTR 


0.5 + RSTRH 


3 + RSTRH 


(172) 


Moment of inertia: 


I = t HSTR [ 2e 2 + - S ^-f - + 
xx s 12 


HSTR 


- e 


+ RSTRH (HSTR - e) 


(173) 
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The equations used to size the stiffener thickness , t g , and 
the area of primary longerons, A LU Lg LL » secondary longerons, 
A__ , and cutout longerons, ALC „, are explained in reference 2. 
A summary of the equations and constraints used in subroutine 
FBEND is given here. The constraints and equations used to 
size the stiffener thickness and longeron areas are: 

1. Minimum area: 


a lu , ls , ll 


> ASTRM 


(173) 


2. Thickness to resits forced crippling: During the 

synthesis of minor frames an approximation to the stiffener area 
based on minimum area and strengths requirements was made. The 
diagonal tension angle, a, was determined for forced crippling 
in the ring. The longitudinal stiffener is now sized to resist 
the forced crippling using the geometry of figure 21 and the a 
computed in MINFR. 

The average stress is defined in equation 160 for stringers 
and longeron construction. Sizing to the allowable stress for 
the stiffener proceeds as in equation 169. 

The maximum induced stress is then 


K f 


FST 


'STmax ” tan a[ t ( 


X B + X C J 


FST 
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where 
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STmax 
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1 + (1 - K) 
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BP AN' 
SFRM 


(174) 


(175) 


The maximum load carried by the stiffeners is 


P 

max 


f e _ A 
STmax s 


(176) 
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3. Maximum down bending carried by upper longitudinal 
members: The nominal maximum allowable stresses for the longi- 

tudinal members and covers are established as: 


Longeron: F „ * 0.9 F 

max cy 

Cover : SC max * 0.76 


(177) 

(178) 


If different materials are used 


F max = min F max' SC max E 




(179) 


and 


SC max = min ( SC max' F max W', 


(180) 


Initial stiffener area estimates are 


Stringers 

Longerons 

A lu = ASTRM 

A LU = AFC 

a ls = AFC 

a ls = 0 

A tt = ASTRM 

A tt - AFC 

LIj 


A Tm = ASTRM 


T T / 

cutout longerons are: 


Co , if RTU = 0 

ALC " 


u 

^ ASTRM, if RTU > 0 


CO , if RTL = 0 

ALC p = 



^ ASTRM, if RTL > 0 


(181) 


(182) 


( 183 ) 
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where RTU and RTL are the apparent width of the section out 
from the shell. 

The down bending moment carried by upper cover, side and 
intermediate longitudinal members is: 
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AV 


CU 


SC max ' 
Z max 


r I i 
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BU 
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r F 'i 

max 1 
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( 184 ) 


where Z . is the distance to the upper cover extreme fiber, 
max 

If 2m av less than the required moment, the difference in 

the moments, AM U , is used to size the longitudinal area 


A LU i 


AM Z 
u max 


F 


max 


h 

AL 


( 185 ) 


If stringer construction is used, an additional cutout 
longeron on each side is inserted to carry the bending loads 
which should have been reacted by stringers. The area is 
sized as: 


ALC u 1 


AM. 

'rtu] 
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Si 

BU i 

max 


2F 


max 


3cu 


( 186 ) 


where Z cu is the vertical distance from the neutral axis to 
the extreme fiber of upper sector. 


4 . The down bending moment carried by the lower sector in 
compression is similar to the tensile side, except only a por- 
tion of the cover material is considered to be effective in 
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carrying load. The critical buckling stress for thin walled 
curved sectors is: 



5 / 3 

4/3 
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9 

'^C 2,1 

a; _ 4 . ,1 i ft 
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where 


WPAN 


bl for longerons 
B gTR for stringers 


B„ d = stringer spacing 

bl K 

if F is greater than the SC the center cover is effec- 
ccr v max 

tive. If not, an effective width of the cover is computed. 
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The moment carried by the effective cover is 


K C - 
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where 
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AV 


BL for longerons 

BL-RTL for stringers 


Z . = vertical distance to extreme fiber of lower sector 

min 
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If SC < F , 
max - ccr' 
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The total moment carried by the ccvvir , side stiffeners, and 
secondary longerons is: 
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M _ m i 1 max 
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The difference between the available moment and the applied 
moment, M t , for the lower cover is: 
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AM, 


ext - M 
2 AV 


(192) 


If AM 0 is greater than zero, the longitudinal areas sized are: 
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(194) 


Cutout Longerons 


Combined Loading 


ALC ■ 


AM £ RTL Z min 
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A LL i 


(F ) (BST) (APL) + P„ v 
max max 


(0-9) (Fey) 


(195) 


The procedure is then repeated for the maximum up-bending case. 
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5. The upper longitudinal members and cover are checked 
tor up-bending assuming all members are in compression. 

6. The lower longitudinal members are checked for up- 
bending assuming all members are in tension. 

7. Vertical stiffness requirement (EI^) : The total verti 

cal bending stiffness available (EIVA) is determined from the 
sum of the components. Because the structure is sized for the 
postbuckled configuration, only those members contributing to 
the shell strength are considered. 
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Longitudinal member stiffness. 


EIVA = 


EIVA C + 
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If the available stiffness is less than the required stiffness 
(EIVT) , the :i fference is used to size the following areas: 


AEI = ETVT - EIVA 


(198) 
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r 1 ) ' 


if RTL < 0 


8. Side bending stiffness requirement (EI Z ) s The total 
side bendiny stiffness available (EISA) is determined from the 
sum of the components. 


EISA o - 1 E chr 
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( 201 ) 
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Longeron 


Stringer 



If the available side stiffness is less than the required stiff- 
ness (EISD) , the difference is used to size the following areas: 

AEI « EISD - EIVA (203) 


The members resized are: 
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The unit weight of longitudinal members is determined as follows 
for the two construction types: 
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Longerons 
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where 


W = widths of upper and lower cutouts 

cut u,* 

BSTR = stringer spacing 
N = number of secondary stringers 

N = number of longerons 

Miscellaneous and Secondary Structural Weight 

The estimation of a number of miscellaneous items under the 
heading Basic Structure (see figure 14) are computed from sta- 
tistical methods. The joints, splices and fasteners are com- 
puted as a fraction of cover and longitudinal member weight. 
Longitudinal partition weight is estimated as a fraction of 
cover and minor frame weight. Flooring and supports items are 
estimated on a unit weight basis for the type of floor used and 
its width and surface area. The engine drag beam and the 
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fitting weights are computed empirically to estimate the attach- 
ment weight penalty for a buried-engine concept and for the 
attachments for the wing, horizontal tail, vertical tail, na- 
celles, and other fuselage mounted components. 

Secondary structural items for which statistical weights 
are computed are generally separated into secondary structural 
items such as canopies, windshields, window, and radomes, and 
into doors, panels, and other access structure. Three methods 
of weight prediction are available for most of these items: a 

rule-of-thumb estimate, a statistical calculation based on 
additional component definition, and direct input by the user. 
The program is arranged on input into three variable groups 
according to calculation requirements. The weight indicator and 
the component CG are required to decide whether to calculate the 
item and its location. If a calculative technique is used, a 
second group of additional required variables must be input. 

The third group of variables may optionally be input in lieu of 
variables in the second group when available. 

The equations used to estimate the secondary structural 
components are defined in reference 2. At least one coefficient 
in each of the statistical equations way also be redefined by 
the user on input. These coefficients are stored in the program 
array, EQN, and are initialized by the program input routine. 
Subsequent input variables may be used to override the default 
values . 


Summary and Intermediate Output 

A summary output of both the computed results and interme- 
diate program values is available. Six print variables currently 
control the output of a'*l print with the exception of certain 
error messages. These suppress or print the input and output 
summaries and control the printing of the intermediate calcula- 
tions from geometry descriptions, and the synthesis of the 
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covers, the bending material, and the minor rings. The inter- 
mediate print is available for most parameters ubed in the 
synthesis loop at each out. These consist primarily of tempo- 
rary variables that are not saved for the final summary print. 

The summary print contains the basic geometric descriptors, 
material gages, and weights computed during the fuselage synthe- 
sis. The basic shell cut geometry and segment properties used 
in the structural design of the bending and torsional stiffnesses 
are given. The shell material gages required to satisfy panel 
flutter, torsion, minimum gage, and strength requirements in 
the skins and landings as well as the required stiffener and 
longeron areas including the effects of cutouts are output for 
each cut. The estimated weight for each segment for each of 
the skin, frames, and longitudinal members is given. A basic 
AN-9102-D weight statement and accompanying balance statement 
as given in figure 14 is used to summarize the weight prediction. 

Optional output also includes a print of the inputs gener- 
ated for the PASS program and plots of the basic configuration 
and material distribution. 

Fuselage Module Sample Case 

A sample test case of the C-141A transport aircraft was run 
as check of the program's accuracy. The computer run was for a 
metallic design with stringer construction. A fixed spacing on 
stringers and minor frames was used, and only the optimum search 
for bulkhead weight was exercised. The input data were derived 
from the demonstration case in reference 2. Since no external 
loads calculation exists within the program, a number of shear 
and bending moment distributions were excerpted from the demon- 
stration case output. The detailed weight statement and balance 
summary are shown in figure 14. 
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The ratio of the C.-141A fuselage weight computed to the 
weight computed here in reference 2 is 1.05. The ratio of the 
actual T141A fuselage to the weight computed in figure 14 is 
1.15. The computed weights include the weight index factors for 
the cover , longerons , frames and bulkheads suggested in refer- 
ence 2 to estimate non-optimum weight increments. About eighty 
percent of the individual components estimated agreed exactly 
with the computations in reference 2. Since major portions of 
the code were adapted directly , the close correlation was expected. 

The major source of differences in estimated structural 
weight was the lack of sufficient loading input. A total of 
sixteen loading conditions were used in the demonstration pro- 
blem. Of these only two occurred in the summary output in 
complete enough form to be used as input in this check case. In 
general , the new geometry descriptors were within a percent of 
the rounded-rectangle representation. The other source of dif- 
ferences was the estimation of the major frames. An empirical 
relation was used in lieu of detailed frame calculations. Though 
the individual frame weights varied significantly, the net 
weight of frames and bulkheads agreed within 3 percent. 

A detailed comparison of each line item in the AN-9102-D 
weight statement was not possible due to the unavailability of 
the C-141A data. A comparison of data from AN-9103-D data is 
given in the following table. The weights are broken down ac- 
cording to the summaries of the items on each of the three pages 
in the detailed weight statement. The comparison of the C-141A 
group data with the results computed by the present program and 
those presented in reference 2 are given in Table IV. 

As group items the estimation of secondary structural com- 
ponents does appear to be satisfactory. The estimation of the 
basic structure still requires additional correlation to estab- 
lish weight index factors, and the proper generation of the 
load spectrum. 


TABLE IV 


Comparison o£ Estimated and Actual C-141A 
Fuselage Group Weight Statements 


Item / Source 

C-141A 

SWEEP 

Fuselage 
Code 

Body Group (total) 

29,342 

26,679 

25,462 

Fuselage or Hull-Basic Structure 

21,438 

18,859 

17,909 

Secondary Structure-Fuselage of Hull 

1,060 

1,138 

1,049 

-Doors, Panels & Misc. 

6,844 

6,680 

6,503 


Conclusions and Recommendations 

The objectives of this work were the improvement of the 
deficiencies in the WADES program for wing design identified in 
reference 1 and the development of a program for the estimation 
of fuselage weights. The enhancements were incorporated into 
the wing design code and comparisons of the results for the F-5A 
were made. A program for the estimation of aircraft fuselage 
weights was written which adapted the descriptors used in the 
ACSYNT program and the methods employed in reference 2. A check 
case for the C-141 was executed and the estimated weights compared. 

The wing design code (shown in both the estimation of the 
material and the load distributions improvement due to the in- 
corporated enhancements) . The extended geometric representation 
of the wing carry-through structural loads and the inclusion of 
trim significantly improved the load distributions. The inclu- 
sion of the empirical calculation of secondary structural com- 
ponent weights and the extended planform definition lowered the 
ratio of computed to structural weight from 1.9 to 1.2 for the 
F-5A. The distributions of structural area and inertia improved 
slightly. The best comparisons were with the surface fit of t, 
and the redesign of only t(E,n). The closest weight estimation, 
however, occurred when redesigning with spar caps. 
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The strength design procedure proved to be sensitive to 
the influence of spar or stringer spacing, bus. che use of the 
free coefficients associated with the sizing of the stringers in 
the synthesis process was inconclusive without the additional 
constraints on rib spacing and global panel buckling. A more 
general set of constraints including core and rib design would 
be required to compute these properly. The inclusion of the 
approximate flutter constraint was very efficient after the ini- 
tial derivative evaluation. The convergence for the combined 
flutter and strength design appeared limited only by the con- 
vergence of the internal load distribution. 

The generation of a computer code for the estimation of 
aircraft fuselage component weights was accomplished. The es- 
timated weights for the check case of the C-141 were within 
sixteen percent of the actual weight. This weight discrepancy 
is due in part to inadequate definition of the load profile. 

No external load calculations were incorporated in the program 
initially. All shear and moment distributions were read in 
directly in lieu of procedure which would generate loads on a 
full aircraft configuration. The coupling of this program to 
an external aerodynamic loads estimating routine is the recom- 
mended extension for improved usability of the program. The 
program was structural initially so that this additional feature 
could be readily incorporated. 
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Figure 2 .- Structural planfortn descriptions. 
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Figure 5.- Equivalent reactions at wing-fuselage junction. 
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Figure 8.- Geometric description of wing-stiffened 
cover sheet model. 
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Figure 23.- Geometric model used in analysis 
of longitudinal members. 
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